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Abstract 



With the solar and SN87a neutrino observations as proofs of concepts, the kilometer-scale neutrino 
experiment IceCube will scrutinize its data for new particle physics. In this paper we review the 
prospects for the realization of such a program. We begin with a short overview of the detector 
response and discuss the reach of "beam" luminosity. After that we discuss the potential of IceCube 
to probe deviations of neutrino-nucleon cross sections from the Standard Model predictions at 
center-of-mass energies well beyond those accessible in man-made accelerators. Then we review 
the prospects for extremely long-baseline analyses and discuss the sensitivity to measure tiny 
deviations of the flavor mixing angle, expected to be induced by quantum gravity effects. Finally 
we discuss the potential to uncover annihilation of dark matter particles gravitationally trapped 
at the center of the Sun, as well as processes occurring in the early Universe at energies close to 
the Grand Unification scale. 
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I. INTRODUCTION 



Conventional astronomy spans 60 octaves in photon frequency, from 10 4 cm radio- waves 
to 10~ 14 cm 7-rays of GeV energy. This is an amazing expansion of the power of our 
eyes which scan the sky over less than a single octave just above 10~ 5 cm wavelength. In 
recent years, detection and data handling techniques of particle physics have been reborn 
in instrumentation to probe the Universe at new wavelengths, smaller than 10~ 14 cm. In 
addition to the traditional photon signals, neutrinos and very high-energy protons (that 
are only weakly deflected by the magnetic field of our galaxy) have become astronomical 
messengers from the Universe. As exemplified time and again, the development of novel ways 
of looking into the cosmos invariably results in the discovery of unanticipated phenomena. 
For particle physicists, the sexiest astrophysical problem is undoubtedly how Nature manages 
to impart an energy of more than one Joule to a single elementary particle. 

Although cosmic rays were discovered almost a century ago, we do not know how and 
where they are accelerated. This may be the oldest mystery in astronomy and solving it is 
challenging as can be seen by the following argument. It is reasonable to assume that, in 
order to accelerate a proton to energy E in a magnetic field B, the size R of the accelerator 
must encompass the gyro radius of the particle: R > Rgy ro = E/B, i.e. the accelerating 
magnetic field must contain the particle's orbit. By dimensional analysis, this condition 
yields a maximum energy E = TBR. The T-factor has been included to allow for the 
possibility that we may not be at rest in the frame of the cosmic accelerator, resulting in the 
observation of boosted particle energies. Opportunity for particle acceleration to the highest 
energies is limited to dense regions where exceptional gravitational forces create relativistic 
particle flows. All speculations involve collapsed objects and we can therefore replace R by 
the Schwarzschild radius R ~ GM / c 2 to obtain E < TBM. 

Cosmic accelerators are also cosmic beam dumps producing secondary photons and neu- 
trino beams. Particles accelerated near black holes pass through intense radiation fields or 
dense clouds of gas leading to production of secondary photons and neutrinos that accom- 
pany the primary cosmic ray beam. The target material, whether a gas or photons, is likely 
to be sufficiently tenuous so that the primary beam and the photon beam are only partially 
attenuated. 

At this point a reality check is in order. Such a dimensional analysis applies to the 
Fermilab accelerator: 10 kilogauss fields over several kilometers (covered with a repetition 
rate of 10 5 revolutions per second) yield 1 TeV. The argument holds because, with optimized 
design and perfect alignment of magnets, the accelerator reaches efficiencies matching the 
dimensional limit. It is highly questionable that nature can achieve this feat. Theorists can 
imagine acceleration in shocks with an efficiency of perhaps 1 — 10%. 

Given the microgauss magnetic field of our galaxy, no structures are large or massive 
enough to reach the energies of the highest energy cosmic rays. Dimensional analysis there- 
fore limits their sources to extragalactic objects. A common speculation is that they may 
be relatively nearby active galactic nuclei powered by a billion solar mass black holes. With 
kilo-Gauss fields we reach 10 11 GeV. The jets (blazars) emitted by the central black hole 
could reach similar energies in accelerating sub-structures boosted in our direction by a 
T-factor of 10, possible higher. The neutron star or black hole remnant of a collapsing su- 
permassive star could support magnetic fields of 10 12 G, possible larger. Shocks emanating 
from the collapse black hole could be the origin of gamma ray bursts and, possibly, the 
source of the highest energy cosmic rays. 
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The astrophysics problem is so daunting that many believe that cosmic rays are not the 
beam of cosmic accelerators but the decay products of remnants from the early Universe, 
for instance topological defects associated with phase transitions (near 10 16 GeV) in Grand 
Unified Theories (GUT's). A topological defect will suffer a chain decay into GUT particles 
X that subsequently decay to familiar weak bosons, leptons, quarks- or gluon-jets. Cosmic 
rays are the fragmentation products of these jets. 

All in all, where the highest energy cosmic rays are concerned, both the accelerator mecha- 
nism and the particle physics are enigmatic. There is a realistic hope that the oldest problem 
in astronomy will be resolved soon by ambitious experimentation. One such experiment is 
the IceCube neutrino telescope, which is required to be sensitive to the best estimates of 
potential cosmic ray neutrino fluxes. Though this telescope is primarily motivated by these 
astronomical goals, it has also appeared in the U.S. Roadmap to Particle Physics, and as 
we will argue in this review, deservedly so. 

As the lightest of fermions and the most weakly interacting of particles, neutrinos occupy 
a fragile corner of the Standard Model and one can realistically hope that they will reveal 
the first and most dramatic signatures of new physics. IceCube's opportunities for particle 
physics are only limited by imagination; they include 

• The search for theories where particle interactions, including gravity, unify at the TeV 
scale. Neutrinos with energies approaching this scale will interact gravitationally with 
large cross sections, similar to those of quarks and leptons, and this increase should 
yield dramatic signatures in a neutrino telescope including, possibly, the production 
of black holes. 

• The search for deviations from the neutrino's established oscillatory behavior that 
result from non-standard interactions, for instance neutrino decay or quantum deco- 
herence. 

• The search for a breakdown of the equivalence principle as a result of non-universal in- 
teractions with the gravitational field of neutrinos with different flavors. Similarly, the 
search for breakdown of Lorentz invariance resulting from different limiting velocities 
of neutrinos of different flavors. 

• The search for neutrinos from the annihilation of dark matter particles gravitationally 
trapped at the center of the Sun. 

• The search from particle emission from cosmic strings or other topological defects and 
heavy cosmological remnants created in the early Universe. 

• The search for magnetic monopoles. 

Alternatively, it is possible that we may be guessing the future while holding too small deck 
of cards and IceCube will open a new world that we did not anticipate. 

The case for doing particle physics with cosmic neutrinos is compelling, the challenge has 
been to deliver the technology to build the instrumentation for a neutrino detector with the 
largest possible effective telescope area to overcome the small neutrino cross section with 
matter, and the best possible angular and energy resolution to address the wide diversity of 
possible signals. We discuss this next. 
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II. OVERALL DETECTOR PERFORMANCE 



In deep ice neutrinos are detected by observation of the Cerenkov light emitted by charged 
particles produced in charged current (CC) and neutral current (NC) interactions. To a 
first approximation, a neutrino (of energy E v ) incident on a side of area L 2 will be detected 
provided it interacts within the lattice of photomultiplier tubes (PMT's) constituting the 
sensitive volume ~ L 3 of the detector. That probability is 

P(E V ) = 1 - exp[-L/l u (E u )] ~ L/l u (E u ) , (1) 

where l u (E v ) = [p icc Na ct ii n{.E ii )\~ 1 is the mean free path. Here p ice = 0.9 g cm -3 is the 
density of the ice, Na = 6.022 x 10 23 is Avogadro's number and a u ^(E u ) is the neutrino 
nucleon cross section. A neutrino flux dF/dE u (neutrinos per GeV per cm 2 per s) crossing 
a detector with energy threshold E^ 1 and cross sectional area A (~ L 2 ) facing the incident 
beam will produce 

f dF 

M = T A{E V ) P{E U ) — dE v (2) 
J £?th arj v 

events after a time T. In practice, the "effective" detector area A is not strictly equal to the 
geometric cross section of the instrumented volume facing the incoming neutrino because 
even neutrinos interacting outside the instrumented volume, may produce a sufficient amount 
of light inside the detector to be detected. Therefore, A is determined as a function of the 
incident neutrino direction by simulation of the full detector, including the trigger. To be 
realistic such a simulation must account for all of the detector properties that we discuss 
next. 

The Antarctic Muon And Neutrino Detector Array (AMANDA) |l] is located below the 
surface of the Antarctic ice sheet at the geographic South pole. During 1993 and 1994, in an 
exploratory phase project, the four- string AMANDA- A array was deployed and instrumented 
with 80 PMT's spaced at 10 m intervals from 810 to 1000 m. A deeper array of 10 strings, 
referred to as AMANDA-B10, was deployed during the austral summers between 1995 and 
1997, to depths between 1500 and 2000 m. The instrument volume of AMANDA-B10 forms 
a cylinder with diameter 120 m, overlooked by 302 PMT's. During December 1997 and 
January 2000, the detector was expanded by adding nine strings of PMT's. The composite 
array of 19 strings and 677 PMT's forms the AMANDA-II array. 

Overall, AMANDA represents a proof of concept for the kilometer-scale neutrino obser- 
vatory, IceCube 0, now under construction. IceCube will consist of 80 kilometer-length 
strings, each instrumented with 60 10-inch photomultipliers spaced by 17 m. The deepest 
module is 2.4 km below the surface. The strings are arranged at the apexes of equilateral 
triangles 125 m on a side. The instrumented (not effective!) detector volume is a cubic 
kilometer. A surface air shower detector, Ice Top, consisting of 160 Auger-style Cerenkov 
detectors deployed over 1 km 2 above IceCube, augments the deep-ice component by provid- 
ing a tool for calibration, background rejection and air-shower physics. The expected energy 
resolution is ±0.1 on a logio scale. Construction of the detector started in the Austral sum- 
mer of 2004/2005 and will continue for 6 years, possibly less. At the time of writing, data 
collection by the first string has begun. 

The event signatures are grouped as tracks, showers, or a combination of the two. Tracks 
include muons resulting from both cosmic ray showers and from CC interaction of muon 
neutrinos. Tracks can also be produced by r leptons arising in ultra-high energy u T CC 
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interactions. Showers are generated by neutrino collisions [y e or low energy v T CC inter- 
actions, and all NC interactions) inside or near the detector, and by muon bremsstrahlung 
radiation near the detector. 



A. Muon Tracks 

In a CC event a produces a muon traveling in nearly the same direction as the neutrino. 
Secondary muons range out over kilometers at ~ 10 3 GeV, to tens of kilometers at 
~ 10 9 GeV, generating showers along their track by bremsstrahlung, pair production 
and photonuclear interactions. All of these are sources of blue Cerenkov light. As the 
energy of the muon degrades along its track, the energy of the secondary showers diminishes 
and the distance from the track over which the associated Cerenkov light can trigger a 
PMT becomes smaller. The geometry of the lightpool surrounding the muon track over 
which single photo-electron are produced is therefore a kilometer-long cone with gradually 
decreasing radius. The orientation of the Cerenkov cone reveals the neutrino direction, 
with an angular resolution »s 0.7° 2]. Muons created by cosmic ray interactions in the 
atmosphere constitute the main background. Full event reconstruction, with quality cuts 
applied for the rejection of cosmic ray muon background, is referred to as "level 2" 0. Down- 
going atmospheric neutrinos at level 2 are highly suppressed because low energy muons need 
to be rejected with energy cuts. 

High energy muons lose energy catastrophically according to 

^ = -a - (3,E, , (3) 

where a = 2.0 x 10~ 6 TeVcm 2 g _1 (characterizes the ionization process) and (3^ = 
4.2 x 10~ 6 cm 2 g~ x (takes into account bremsstrahlung, e + e~ pair production and nuclear 
interactions) 0. The distance a muon travels before its energy drops below some energy 
threshold _E* h , called the muon range, is then given by 



Z min = 1 ln 



11 P 



(4) 



In the first kilometer a high energy muon typically loses energy in a couple of showers of 
one tenth its energy. So the initial size of the cone is the radius of a shower with 10% of the 
muon energy, e.g. 130 m for a 100 TeV muon. Near the end of its range the muon becomes 
minimum ionizing emitting light that creates single photoelectron signals at a distance of 
just over 10 m from the track. For 0.3 photoelectrons, the standard PMT threshold setting, 
this distance reaches 45 m. 

The expected event rate at IceCube can be estimated through semi-analytical calcu- 
lations 

pi poo poo 

fin / J ttiO 



M Vfl = T / dcos9 / dip dEf \ dE° / dE v (5) 

'"' " ! ' (E,, cos 9)^(E v ,E")rvrF (El E^,QA° eS , 



dE v d cos9" ' ' dE° 



where dF u ^/dE v d cos 6 is the differential muon neutrino neutrino flux, dacc/dE® is the 
differential CC interaction cross section producing a muon of energy E® and ut is the 
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FIG. 1: Effective function of the final muon energyfin GeV) after level 2 cuts. The full 

line indicates the result from the semi-analytical calculation [5|], whereas the data points are from 
Monte Carlo simulation [fj. For comparison, also shown by a dashed line is Ao(E^ n ). 




FIG. 2: Expected M v event rate after level 2 cuts for one year exposure obtained with the 
semianalytical calculation (full lines) [jj and with Monte Carlo simualtions (dashed lines) Q]. 
The left panel shows the energy spectrum (E v in GeV) and the right panel the zenith angular 
distribution. 



number density of nucleons in the matter surrounding the detector. After production with 
energy E®, the muon ranges out in the rock and in the ice surrounding the detector and 
looses energy. Here, F(E®, E^ n , / M ) is the function that describes the energy spectrum of 
the muons arriving at the detector. Hence, F(E®,E^ n ,ln) represents the probability that 
a muon produced with energy E® arrives at the detector with energy E^ n after traveling 
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a distance l^. The function F(E®, E^ n , / M ) is computed by propagating the muons to the 
detector taking into account energy losses as given in Eq. (jlj). The possibility of fluctuations 
around the average muon energy loss (using the average energy loss would equalize to 
the average muon range distance) is included in F(E®, Thus, E®, E^ n , and, Z M are 

kept as independent variables. For simplicity, and F(E® E^ n , l^j in ice are used and the 
effect of the rock bed below the ice is accounted for in the form of an additional angular 
dependence of the effective area for upward going events. The details of the detector response 
are encoded in the effective area A° cE . A phenomeno logical parametrization has been used 
to simulate the response of the IceCube detector after events that are not neutrinos have 
been rejected (i.e., level 2 cuts [6]) 



A^ = ME??) x ^( cos ^ E T) x R ( l 7 n ) > ( 6 ) 

where A G (E^) includes the energy dependence of the effective area due to trigger re- 
quirements 6]. Here i?(/™ m ) represents the smearing in the minimum track length cut, 
I™ = 300 m, due to the uncertainty in the track length reconstruction which can be 
parametrized by a Gaussian 

i poo (lmin ^/min\2 

mf n ) = ^^ dff^exp- ^ J } , (7) 

\l L-KOl Jo Zcr l 

with <Ji = 50 m. Introduction of a simple straight line dependence on log 10 (.E* n ) 



A (E^) = Ao 



1 + 0.55 log 



10 



GeV 



(8) 



leads to a good agreement with the results of Monte Carlo simulations |6j. Here Ao is an 
overall normalization constant, which is calibrated using the flux of atmospheric neutrinos 
as seen by AMANDA 0. After level 2 cuts are applied the atmospheric neutrino flux yields 
91000 events/yr at IceCube. Next, one has to "simulate" the cuts in the muon tracklength 
/™ m and the number of PMT's reporting signals in an event p. The angular dependence 
of the effective area for downgoing events (8 < 80°) is determined by the level 2 cut on the 
minimum number of channels Ncn > N^(cos6) = 150 + 250 cos 9. This requirement leads 
to an ^"-dependent angular constraint 



R(cos o, eT) = * r dNcH exp j^ch-(w _ (9) 

V2Trcr NcH J N ™ Z(T n ch 

where (Ncn) E 6n is the average channel multiplicity produced by a muon which reaches the 
detector with energy E? n and <J% CK is the spread on the distribution. The parametrization 

X 



log 10 (Vch)^) = 2.0 + 0.88-== , (10) 



with 

X = 0.47 (log u (|£ )-4.6 | (11) 
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and 

a NcH =OA(N CR ) E ^ (12) 

shows good agreement with the charged multiplicity obtained through Monte Carlo sim- 
ulations p. Finally, to account for the presence of the rock bed below the detector, one 
can introduce a phenomenological angular dependence of the effective area for upward going 
muons 

^(cos 9) = 0.70 - 0.48 cos 9 for 9 > 85° , (13) 

independent of the muon energy. Figure d shows the effective area A e s, defined as the ratio 
of the number of upgoing muon events, with/without the inclusion of Ao(E^ n ) x R(l™ m ) and 
the level 2 cuts on Z™ m , and compare the results of the semi-analytical calculation with Monte 
Carlo simulations of the detector response . Figure 121 compares the energy spectrum and 
the zenith angular distribution of the events (after level 2 cuts are applied) obtained with 
the semi-analytical calculation and with Monte Carlo simulations. 



B. Electromagnetic Showers 

Showers are generated by neutrino collisions — v e or v T CC interactions, and all NC 
interactions — inside of or near the detector, and by muon bremsstrahlung radiation near 
the detector. Normally, a reduction of the muon bremsstrahlung background is effected by 
placing a cut of 4 x 10 4 GeV on the minimum reconstructed energy j8j. 

Electron neutrinos deposit 0.5-0.8% of their energy into an electromagnetic shower ini- 
tiated by the leading final state electron. The rest of the energy goes into the fragments 
of the target that produce a second subdominant shower. Because the size of the shower, 
of order meters in ice, is small compared to the spacing of the PMT's, it represents, to a 
good approximation, a point source of Cerenkov photons radiated by the shower particles. 
These trigger the PMT at the single photoelectron level over a spherical volume whose radius 
scales linearly with the shower energy. For ice, the radius is 130 m at 10 4 GeV and 460 m 
at 10 10 GeV, i.e. the shower radius grows by just over 50 m per decade in energy. 

The measurement of the radius of the sphere in the lattice of PMT determines the en- 
ergy and turns neutrino-detection-experiments into total energy calorimeters. Note that a 
contained "direct hit" by a 10 10 GeV neutrino will not saturate a km 3 detector volume. 
Because the shower and its accompanying Cerenkov lightpool are not totally symmetric but 
elongated in the direction of the leading electron (and incident neutrino), its direction can be 
reconstructed. Pointing is however far inferior to what can be achieved for muon neutrinos 
and estimated to be precise to ~10° only. 

The expected shower event rate at IceCube can also be estimated through semi-analytical 
calculations 0, 

Mh = A/" S h,cc + A/" sh ,Nc , (14) 

where 

A/"sh,cc 

and 



Tn T / dE v V -—(E v )acc(Eu)Vos(E v ) , (15) 



poo poo rl~F l ' a drr 

Mh,Nc = Tn T / dE'J dE„ -JW- m-^i^EttVrtiK) ■ (16) 

J E,,—E n ? [n JE™ in * at ^v atlj v 
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FIG. 3: Effective volume of IceCube for showers initiated by neutrinos. The parametrization shown 
by a solid curve corresponds to V eS {E u ) « -0.071ogf E„ + 1.271og 10 ^ - 3.21. 

Here, d<T NC /dE' u is the differential NC interaction cross section producing a secondary neu- 
trino of energy E' v , and V e a(£v) is the effective volume shown in Fig. 01 In writing Eqs. fl"5|) 
and we are assuming that for contained events the shower energy corresponds with the 
interacting u e or v T neutrino energy (E^ = E u ) in a CC interaction, while for NC the shower 
energy corresponds to the energy in the hadronic shower E U = E v — E' v = E v y, where y is 
the usual inelasticity parameter in deep inelastic scattering |l0j . 




C. Shower- Track Combo 

For z/ r 's, CC current interactions have different signals depending on the energy. The r 
has a decay length 

I T ^50mx {E T /10 6 GeV) , (17) 
and loses energy in the Earth according to 

^ = -a-(3 T E T , (18) 

where (3 T = 8 x 10~ 7 cm 2 g _1 0]. From Eqs. (j!7)l and ()18j) it is easily seen that for energies 
above about 10 8 GeV the r is more likely to interact electromagnetically than to decay, 
whereas for E T < 10 7,5 GeV the energy loss before decay is negligible. 

For t leptons less energetic than 10 6 GeV, the shower (hadronic or electromagnetic) from 
the tau decay cannot be separated from the hadronic shower of the initial v T interaction. 
At a few times 10 6 GeV, the range of the tau becomes a few hundred meters and the two 
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showers produced may be easily separated and be identify as a double bang event. At 
energies between 10 7 — 10 7,5 GeV, the tau decay length is comparable to the instrumented 
volume. In such cases, one may observe a r track followed by the r-decay shower ("lollipop 
topology"), or a hadronic shower followed by a r track which leaves the detector ("popillol 
topology"). At energies > 10 7 ' 5 GeV, l T 3> 1 km and r's leave only a track like muons. 
However, a r going through the detector at high energies without decaying will not deposit 
as much energy in the detector as a comparable-energy muon, due to the mass difference. 
(The direct-pair production process scales inversely with mass, so it dominates tau-lepton 
energy loss resulting in l/20 th the light produced by a muon.) Such a tau might then 
be indistinguishable from a low energy muon track, and thus might not be flagged as an 
interesting event. In summary, the energy range from 10 6 — 10 7 ' 5 GeV is the "sweet spot" 
for r detection in IceCube, since here one can observe all the distinctive topologies. 



Although neutrino "telescopes" are designed as discovery instruments, be it for particle 
or astrophysics, their conceptual design is very much anchored to the observational fact that 
Nature produces protons and photons with energies in excess of 10 11 GeV and 10 4 GeV, 
respectively. The cosmic ray connection sets the scale of cosmic neutrino fluxes. We discuss 
this next. 



III. LUMINOSITY 

IceCube is unique among high energy neutrino experiments in being sensitive to ten 
decades in energy. Unlike terrestrial experiments with a man-made beam, the luminosity 
varies significantly with the center-of-mass energy In what follows we discuss the 

different energy regimes accesible to IceCube. 



A. Atmospheric Neutrinos 

A guaranteed beam of high energy neutrinos originates in the atmospheric cascades initi- 
ated by cosmic rays 0|. When protons and nuclei enter the atmosphere, they collide with 
the air molecules and produce all kind of secondary particles, which in turn interact or decay 
or propagate to the ground, depending on their intrinsic propeties and energies. In the GeV 
range, the most abundant particles are neutrinos. The production mechanism is the decay 
chains of mesons created in these cascades. 

Pion decay dominates the atmospheric neutrino production: tt + — > — > e + v e v l j? ll 
and the conjugate process. This decay chain determines the neutrino energy spectra up to 
about 100 GeV above which they become increasingly modified by the kaon contribution, 
which asymptotically reaches 90%. In the atmosphere mesons encounter the interaction- 
decay competition. Therefore, neutrinos from meson decay will have a spectrum one power 
of energy steeper than the primary cosmic ray spectrum. The muon daughter neutrinos 
will have a spectrum steeper by two powers of energy, because the muon spectrum itself is 
steeper by 1/E. Electron neutrinos have a differential spectrum (approximately) oc E~ 4 ' 7 . 
The muon neutrino spectrum is flatter, 

, p u, , , p .„„ _ / 1.6 x lO- 16 ^ 3 - 7 GeVcm-^sr" 1 E v < 10 5 GeV , , 

af atm /au u - <p atm ~ < 5 ^ Q x 1Q _ 18 4.o GeVcm - 2 s -i gr -i E > 1Q 5 GeV • UJJ 



11 



In this energy window, the flavor ratios are w e : : w T « 1/20 : 19/20 : [14| and, as can 
be seen in Fig. HJ the energy spectra are a function of the zenith angle of the atmospheric 
cascades. This is because mesons in inclined showers spend more time in tenuous atmosphere 
where they are more likely to decay rather than interact. For this reason the spectra of highly 
inclined neutrinos are flatter than those of almost vertical neutrinos. In the sub-GeV energy 
range the spectra are significantly flatter (parallel to the primary cosmic ray spectrum) as 
all mesons and muons decay. Thus, one can immediately predict that the flavor ratios at 
production are 1/3 : 2/3 : 0. Above about 10 5 GeV, kaons are also significantly attenuated 
before decaying and the "prompt" component, arising mainly from very short-live charmed 
mesons (D ± , D°, D s and A c ) dominates the spectrum Such a "prompt" neutrino flux 
is isotropic with flavor ratios 12/25 : 12/25 : 1/25. 
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FIG. 4: Energy spectrum of atmospheric neutrinos as seen by the AMANDA II and Frejus [lfi | 
experiments. The two dotted curves indicates model predictions 01 for the horizontal (upper) 
and vertical (lower) flux. 

The neutrino flux arising from pion and kaon decay is reasonable well understood, with an 
uncertainty in the range 10% — 20% The prompt atmospheric neutrino flux, however, is 



much less understood, because of uncertainty about cosmic ray composition and relatively 
poor knowledge of small-x QCD processes [l9j]. For numerical estimates, throughout this 
review we have chosen to use the 3-dimensional estimates of conventional atmospheric neu- 
trino flux given in |2^, which we extrapolate to match at higher energies the 1-dimensional 
calculations of Volkova • We also incorporate "prompt" neutrinos from charm decay as 
calculated in [2l| . 

In 10 yr of operation IceCube will collect more than 7 x 10 5 atmospheric neutrino events 
with El n > 100 GeV. These events are generated by neutrinos with large enough energy 
for standard mass-induced oscillations to be very much suppressed so they should behave as 
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flavor eigenstates. As we discuss in Sec. IV A[ this high-statistics high-energy event sample 
offers a unique opportunity to probe new physics mechanisms for leptonic flavor mixing. 



B. Extraterrestrial Neutrinos 

Interactions of ultra-high energy cosmic ray protons propagating over cosmological dis- 
tances with the cosmic microwave background (CMB) generates a nearly guaranteed cosmo- 
genic flux of neutrinos 0] through the decay of charged pions produced in jry interactions, 
which should also result in a suppression of the cosmic ray spectrum above the "GZK cutoff" : 
Egzk ~ 10 10 ' 7 GeV [23]. The intermediate state of the reaction P7cmb —* Ntt is dominated 
by the A + resonance, because the n decay length is smaller than the nucleon mean free path 
on the relic photons. Hence, there is roughly an equal number of ir + and 7r°. Gamma rays, 
produced via 7r° decay, subsequently cascade electromagnetically on the cosmic radiation 
fields through e + e~ pair production followed by inverse Compton scattering. The net result 
is a pile up of 7 rays at GeV energies, just below the threshold for further pair production. 
On the other hand, each 7r + decays to 3 neutrinos and a positron. The e + readily loses its 
energy through synchrotron radiation in the cosmic magnetic fields. The neutrinos carry 
away about 3/4 of the ir + energy, and therefore the energy in cosmogenic neutrinos is about 
3/4 of the one produced in 7-rays. 

The normalization of the neutrino flux depends critically on the cosmological evolution of 
the cosmic ray sources and on their proton injection spectra 24[ 25|. It also depends on the 



assumed spatial distribution of sources; for example, relatively local objects, such as sources 
in the Virgo cluster |26j|, would dominate the high energy tail of the neutrino spectrum. 
Another source of uncertainty in the cosmogenic neutrino flux is the energy at which there 
is a transition from Galactic to extragalactic cosmic rays as inferred from a change in the 
spectral slope (22|. A fourth source of uncertainty in the cosmogenic flux is the chemical 
composition — if ultra-high energy cosmic rays are heavy nuclei rather than protons the 
corresponding cosmogenic neutrino flux may be somewhat reduced |28| . 

In addition to being produced in the propagation of ultra-high energy cosmic rays, neu- 
trinos are also expected to be generated in their sources. How many neutrinos are produced 
in association with the cosmic ray beam? The answer to this question, among many others, 
provides the rationale for building kilometer-scale neutrino detectors. We first consider a 
neutrino beam produced at an accelerator laboratory. The accelerator beam is "dumped" 
into a large target, for instance a hundred meters of steel. The target absorbs all parent 
protons as well as the secondary electromagnetic and hadronic showers. Only neutrinos, 
the decay products of charged pions, exit the dump. If Nature constructed such a "hidden 
source" in the heavens, conventional astronomy will not reveal it. It cannot be the source 
of the cosmic rays, however, because in this case the dump must be transparent to protons. 

A more generic "transparent" source can be imagined as follows: protons are accelerated 
in a region of high magnetic fields near a black hole or neutron star. They inevitably interact 
with the radiation surrounding collapsed objects via the processes p + 7 — > A — > 7T° + p 
and p + 7-^A^7r + + n. While the secondary protons may remain trapped in the 
acceleration region, equal numbers of neutrons, neutral and charged pions escape/decay. 
The energies escaping the source in the form of cosmic rays, gamma rays and neutrinos 
produced by the decay of neutrons and neutral and charged pions, are related by the physics 
of photoproduction. The neutrino flux from a gen eric transparent cosmic ray source is often 
referred to as the Waxman-Bahcall (WB) flux 29]. It is easy to calculate and the derivation 
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is revealing. 

We will concentrate on the population of extragalactic cosmic rays with energies in excess 
of 10 99 GeV, i.e., beyond the "ankle" in the spectrum. The flux above the ankle is often 
summarized as "one 10 10 GeV particle per kilometer square per year per steradian". This 
can be translated into an energy flux 

E (dN CR } 10 10 GeV 



dE J (10 10 cm 2 )(3 x 10 7 s)sr 

= 3 x 1(T 8 GeV cm" 2 s" 1 sr" 1 . (20) 

From this we can derive the energy density pe in cosmic rays using flux = velocity x density, 
or 

JN C r 



4tt / dE\E-^\=cp E . (21) 



We obtain 

4tt 3 x 1(T 8 _ GeV _ iq TeV 

— ~ 10 

l 3 r»m 3 



Pe = — I "V dE ^ ~ 10- 19 ^ , (22) 
c J F . E cur cur 

Kf J_/ n 



taking the extreme energies of the accelerator(s) to be E max / E min ~ 10 3 . 

This energy density derived more carefully by integrating the spectrum beyond the 
"ankle" assuming an E~ 2 energy spectrum with a GZK cutoff is closer to ~ 3 x 
10 -19 erg cm~ 3 30] . The power required for a population of sources to generate this en- 
ergy density over the Hubble time of 10 10 years is ~ 3 x 10 37 erg s _1 per (Mpc) 3 or, as often 
quoted in the literature, ~ 5 x 10 44 TeV per year per (Mpc) 3 . This works out to |3l| 

• ~ 3 x 10 39 erg s -1 per galaxy, 

• ~ 3 x 10 42 erg s -1 per cluster of galaxies, 

• ~ 2 x 10 44 erg s -1 per active galaxy, or 

• ~ 2 x 10 52 erg per cosmological gamma ray burst. 

The coincidence between these numbers and the observed output in electromagnetic energy 
of these sources explains why they have emerged as the leading candidates for the cosmic 
ray accelerators. The coincidence is consistent with the relationship between cosmic rays 
and photons built into the "transparent" source. In the photoproduction processes roughly 
equal energy goes into the secondary neutrons, neutral and charged pions whose energy 
ends up in cosmic rays, gamma rays and neutrinos, respectively. For E max = 10 12 GeV, 
the generic source of the highest energy cosmic rays produces a flux of E v 2 dN v / dE u ~ 
4 x 10 -8 GeVcm" 2 s _1 sr _1 , close to the back-of-the-envelope calculation performed earlier. 
There are several ways to modify this simple prediction: 

• The derivation fails to take into account the fact that there are more cosmic rays in 
the Universe producing neutrinos than observed at Earth because of the GZK-effect 
and neglects evolution of the sources with redshift. This increases the neutrino flux 
by a factor dn/dcMB, the ratio of the Hubble to the attenuation length radius over 
which cosmic rays above the ankle penetrate the cosmic microwave background. 
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• For proton-7 interactions muon neutrinos (and muon antineutrinos receive only 1/2 of 
the energy of the charged pion in the decay chain tt + — > ji + + — > e + + v e + + 
assuming that the energy is equally shared between the 4 leptons. Furthermore half 
the muon neutrinos oscillate into tau neutrinos over cosmic distances. 

In summary, 

dN v 1 1 dN C R d H dNcK , OQ , 

E v = - x - x E x ~ E (23) 

dE v 2 2 dE c/cmb dE y ' 

The corrections approximately cancel. The transition from galactic to extragalactic sources 
is debated; a transition at lower energy significantly increases the energy in the extragalactic 
component. This raises the possibility of an increase in the associated neutrino flux (32^ . 
Waxman and Bahcall referred to this as a bound because in reality more energy is transferred 
to the neutron than to the charged pion in the source; in the case of the photoproduction 
reaction p + 7 — > A — > 7r + + n it is four times more. Therefore 

dN u _1 dN CR 
Ev dE v ~ \ E ^E ■ (24) 

This is, for instance, the diffuse flux expected from all gamma ray bursts. 

In the end we estimate that for each neutrino flavor the flux associated with the sources 
of the highest energy cosmic rays is loosely confined to the range, 

#^wb ^ 1 .3 x 1(T 8 GeV cm" 2 s" 1 sr" 1 , (25) 

yielding ~ 10 — 50 detected muon neutrinos per km 2 per year. This number depends weakly 
on E max and the spectral slope. The observed event rate is obtained by folding the predicted 
flux with the probability that the neutrino is actually detected in a high energy neutrino 
telescope; only one in a million is at TeV energy. As we discussed in Sec. Ill A[ this probability 
is given by the ratio of the muon and neutrino interaction lengths in the detector medium, 

This flux has to be compared with the sensitivity of ~10~ 7 GeV cm -2 s _1 sr _1 reached 
durin g th e first 4 years of operation of the completed AMANDA detector in 2000- 
2003 |33|. The analysis of the data has not been completed, but a preliminary limit 
of 2.9 x 10~ 7 GeV cm -2 s" 1 sr^ 1 has been obtained with a single year of data [34j. On 
the other hand, after three years of operation IceCube will reach a diffuse flux limit of 
Etch 

iceCube ~ 2 7 x 10 Ge V cm 2 s * sr The exact value depends on the magnitude 
of the dominant high energy atmospheric neutrino background from the prompt decay of 
atmospheric charmed particles p. As noted in the previous section, the level of this back- 
ground is difficult to anticipate. A cosmic flux at the " Waxman-Bahcall" level will result in 
the observation of several hundred neutrinos in IceCube 0. 

In Fig. we plot the expected spectrum of ultra-high energy cosmic neutrinos in the 
models discussed above. For comparison we show the upper limits on neutrino fluxes re- 
ported by the AMANDA Collaboration j3(| and the cascade limit. 1 Hidden sources must 
have a steeper spectrum than transparent sources in order to be consistent with data. In 



The cascade limit arises from the requirement that the diffuse 7-ray fluxes initiated by photons and 
electrons from pion decays should not exceed measurements. The cascade limit from Ref. [i^ exploits 
the measurement of the diffuse 7-ray background by EGRET ; 3S]. A lower extragalactic contribution 
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FIG. 5: Neutrino fluxes from transparent and hidden sources. The horizontal solid line indicates the 
WB prediction (summed over all flavors) which corresponds to a Galactic/extra-galactic crossover 
energy ~ 10 10 GeV [2^. The falling solid line indicates the expected neutrino flux from transparent 
sources, if one assumes the onset of dominance by the extra-galactic component is at 10 8,6 GeV 32]. 
The dash and dashed-dotted lines indicate neutrino flux estimates from hidden sources (HfJ. The 
label on each curve indicates the time required for IceCube to achieve an integrated luminosity 
C ps 30 nb" 1 at y/s ~ 6 TeV. The cross-hatched region excludes an E u 2 spectrum at the 90% 
CL by measurements of AMANDA-BIO [3^ . The single hatched region, obtained by rescaling the 
AMANDA integrated bolometric flux limit to an E~ 2 ^ A power law, is the exclusion region for the 
low crossover model. The shaded region indicates the cascade limit. 



particular, the flux could be enhanced by more than an order of magnitude in the energy 
decade between 10 7 ~ 8 GeV. Therefore, the energy region E u ~ 10 7-7 ' 5 GeV is particularly 
amenable to probe high energy neutrino interactions: commonly proposed non-atmospheric 
neutrino fluxes are still sizeable, while the anisotropic background from atmospheric neu- 
trinos has become insignificant. This insures a statistically significant sample of neutrino 
collisions with center-of-mass energy y/s ~ 6 TeV. It should be noted that this energy is well 
above the HERA domain (y/s ~ 500 GeV), the highest accelerator energy at which even 
indirect tests of the Standard Model cross sections are possible. 



C. Directional Signals 



MeV neutrino astronomy has been possible for about 40 years. Thus far two sources have 
been identified: the Sun |4lJ and supernova 1987a (SN87a) 42]. During the next decade, 



to the 7-ray back gro und than that inferred in above reference, by roughly a factor of two, has been 
proposed recently [39j. The cascade limit may therefore be stronger by a corresponding factor. Note that 
logarithmic corrections to the spectrum oc E~' 2 could result in sizeable deviations of the cascade limit at 
very high energy [dof . 
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the large sensitivity of IceCube to extraterrestrial neutrino fluxes gives a hope to open a 
new era in neutrino astronomy through detection of high energy neutrinos. 

One of the key features of high energy neutrinos is that they are unique carriers of un- 
ambiguous information about hadronic process in cosmic accelerators 1431 . Since generally 
these processes also result in TeV gamma-rays of comparable fluxes j44|, the best candi- 
dates to be discovered by neutrino telescopes are sources of gamma-rays with fluxes above 
1 TeV around (3 — 10) x 10~ 12 phcm~ 2 s _1 . Presently several TeV gamma-ray sources, in 
particular persistant galactic sources like Crab Nebula jiH, shell type SNRs RXJ1713 [46| 
and Vela Jr 14711 . as well as strongly variable Active Galactic Nuclei Mrk 421, Mkn 501, 
and 1ES 1959 |48J] do show fluxes of TeV gamma-rays above 10 -11 phcm -2 s _1 . Therefore, if 
significant fractions of these gamma-ray fluxes are of hadronic origin IceCube should be able 
to detect the neutrino counterparts of TeV gamma-rays. Moreover, in the case of sources 
with heavily absorbed TeV gamma-ray emission, the chances of detection of countrerpart 
neutrino emission could be quite high even for relatively weak TeV sources with a flux well 
below the level 0.1 Crab 0. 

Antineutrinos are also produced through the /?-decay process, but since cosmic neutrino 
production is dominated by 7r ± decay, the search for the V e channel is not straightforward. A 
potential way is to identify neutron-emitting-sources. However, the decay mean free path of 
a neutron is cT n T n = 9.15 (E n /10 9 GeV) kpc (the lifetime being boosted from its rest-frame 
value, ~T n = 886 seconds, to its lab value by T n = E n /m n ), and so only neutrons with energy 
> 10 9 GeV have a boosted cr n sufficiently large to serve as Galactic messengers. 

Air shower arrays have observed a "directional" flux of cosmic rays from the galactic 
plane jH3], unlikely to be protons since their directions are scrambled in the Galactic magnetic 
field. The flux appears only in a narrow energy range from 10 8 ' 9 to 10 9 ' 5 GeV, the energy 
where neutrons reach typical galactic kiloparsec distances within their lifetime of minutes. 
Both the directionality and the characteristic energy make a compelling case for electrically 
neutral neutron primaries. The galactic plane excess, which is r oug hly 4% of the diffuse 
flux, is mostly concentrated in the direction of the Cygnus region [51| . 

Independent evidence may be emerging for a cosmic accelerator in the Cygnus spiral 
arm. The HEGRA experiment has detected an extended TeV 7-ray source, J2032+4130, 
in the Cygnus region with no clear counterpart and a spectrum not easily accommodated 
with synchrotron radiation by electrons [H^- The difficulty to accommodate the spectrum by 
conventional electromagnetic mechanisms has been exacerbated by the failure of CHANDRA 
and VLA to detect X-rays or radiowaves signaling acceleration of any electrons The 
model proposed is that of a proton beam, accelerated by a nearby mini-quasar or possibly 
Cygnus X-3, interacting with a molecular cloud to produce pions that are the source of the 
gamma rays. Especially intriguing is the possible association of this source with Cygnus 
OB2, a cluster of more than 2700 (identified) young, hot stars with a total mass of ~ 10 4 
solar masses j54|. Proton acceleration to explain the TeV photon signal requires only 0.1% 
efficiency for the conversion of the energy in the stellar wind into cosmic ray acceleration. 
Also, the stars in Cygnus OB2 could be the origin of time-correlated, clustered supernova 
remnants forming a source of cosmic ray nuclei. By cooperative acceleration their energies 
may even exceed the ~ 10 6 GeV cutoff of individual remnants and accommodate cosmic rays 
up to the ankle, where the steeply falling (oc E~ 3A6±om ) cosmic ray spectrum flattens to 
^-2.s±o.3 j^| ^ n i mmec li a te consequence of this nucleus- dominance picture is the creation 
of free neutrons via nuclei photodisintegration on background photon fields. These liberated 
neutrons are presumably responsible for the observed directional signals. 
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The observed spectrum from the Cygnus region can be described by single power law 
reflecting the average shape of the diffuse cosmic ray spectrum between 10 6 and 10 8 5 GeV, 



dF n dF n 



e -D/{cV n T n ) 



source 



dE n dE n 

= cr £ r3.i e -B/(cr n r n ) j (26 ) 

where D « 1.7 kpc is the distance to Cygnus OB2. By integrating the spectrum between 
Ei = 10 8 ' 9 GeV and Ei = 10 9 ' 5 GeV jHH, normalization to the observed integrated flux (5ll |. 

/ CE- 3A e- D/{cVnTn) dE n « 9 km~V _1 , (27) 

J E\ 

leads to C = 1.15 x 10 20 km" 2 yr" 1 . 

For every surviving neutron at ~ 10 9 GeV, there are many neutrons at lower energies 
that do decay via n — > p + e~ +V e . The proton is bent by the Galactic magnetic field and the 
electron quickly loses energy via synchrotron radiation, but the V e travels along the initial 
neutron direction, producing a directed TeV energy beam which is potentially observable. 

The basic formula that relates the neutron flux at the source to the antineutrino flux 
observed at Earth is 56l |: 



dF- f dF 

(Ejt) = dE n ^(E n 



dEp y J dE, 
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1 - e B n r n ) I fa- —(e T 



5 [Ep - E n e v (1 + cos 9 v )/m n ] . (28) 

The variables appearing in Eq. (|2*Hj) are the antineutrino and neutron energies in the lab 
{E„, E n ), the antineutrino angle with respect to the direction of the neutron momentum in 
the neutron rest-frame {9jt), and the antineutrino energy in the neutron rest-frame (e^). The 
last three variables are not observed by a laboratory neutrino-detector, and so are integrated 
over. The observable E v is held fixed. The delta-function relates the neutrino energy in the 
lab to the three integration variables. The parameters appearing in Eq. (|2*Hj) are the neutron 
mass and rest-frame lifetime {m n and r ra ). Finally, dPjde v is the normalized probability 
that the decaying neutron produces an antineutrino with energy e v in the neutron rest- 
frame. Note that the maximum antineutrino energy in the neutron rest frame is very nearly 
Q = m n — m p — m e = 0.71 MeV. 

The integral neutrino flux, F F (> Ep) = J E _dEvdFv/dEv, is particularly useful for ex- 
periments having a neutrino detection-efficiency that is independent of neutrino energy, or 
nearly so. The integral flux, normalized to the integrated neutron flux in Eq. (|27|). is shown 
in Fig. |U1 Note that the nuclear photodisintegration threshold implies an infrared cutoff 
on the primary neutron energy at the source, which in turn leads to a low energy cutoff of 
O (TeV) on the integral flux. 



In summary, the AMANDA experiment using natural 1 mile-deep Antarctic ice as a 
Cerenkov detector is in steady operation collecting roughly 7 to 10 neutrinos per day. Ex- 
pansion into the IceCube kilometer scale neutrino telescope will increase the event rate at 
least to 250 events per day. Because of this high-statistics data sample IceCube will be a 
powerful high-energy physics laboratory. We explore the science reach next. 
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FIG. 6: Integrated flux of Vn + V e + v T (solid line) predicted to arrive at Earth from the direction 
of the Cygnus region. Also shown are the integrated + Vn and v e + V e atmospheric fluxes for an 
angular bins of 1° x 1° and 10° x 10°, respectively. The shaded band indicates the region excluded 
by the AMANDA experiment jH^J.The fluxes of neutrinos inferred from HEGRA measurements of 
the 7-ray flux are also shown: the lower line is based on the assumption of pj interactions, whereas 
the upper line is based on pp interactions (the charged/neutral pion-production ratio depends on 
the interaction). In each case the solid portion of the line indicates the region where HEGRA data 
is available and the dashed part is an extrapolation to unobserved energies 



IV. NEUTRINO INTERACTIONS BEYOND THE WEAK SCALE 



The saga of the Standard Model [58| is still exhilarating because it leaves all questions 
of consequence unanswered. The most evident of unanswered questions is why the weak 
interactions are weak. Though electromagnetism and weak interactions are unified, elec- 
tromagnetism is apparent in day life while the weak interactions are not. Already in 1934 
Fermi |59J provided an answer with a theory that prescribed a quantitative relation between 
the fine structure constant and the weak coupling, Gf ~ at/Myy. Although Fermi adjusted 
Myy to accommodate the strength and range of nuclear radioactive decays, one can readily 
obtain a value of Mw of 40 GeV from the observed decay rate of the muon for which the 
proportionality factor is 7r /y/2. The answer is off by a factor of 2 because the discovery of 
parity violation and neutral currents was in the future and introduces an additional factor 
1 - M^/MJ 



1, 



G f 



1-M^/MfJ 



l + Ar). 



(29) 
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Fermi could certainly not have anticipated that we now have a renormalizable gauge theory 
that allows us to calculate the radiative correction Ar to his formula. Besides regular higher 
order diagrams, loops associated with the top quark and the Higgs boson contribute. 

If one calculates the radiative corrections to the mass M# appearing in the Higgs poten- 
tial, the same theory that withstood the onslaught of precision experiments at LEP/SLC 
and the Tevatron yields a result that grows quadratically 

6M h = TTT^ ( 2M w + M 2 Z + M 2 H - AM 2 ) A 2 , (30) 

where M\ = 2\v 2 , A is the quadratic Higgs coupling, v = 246 GeV, and A is a cutoff. Upon 
minimization of the potential, this translates into a dangerous contribution to the Higgs 
vacuum expectation value which destabilizes the electroweak scale j6(|. The Standard Model 
works amazingly well by fixing A at the electroweak scale M\y. It is generally assumed that 
this indicates the existence of new physics beyond the Standard Model. Following Weinberg, 

£(M W ) = ~ M%H*H + I A(H^H) 2 + + £™ awa + - C 5 + + . . . , (31) 

where the operators of higher dimension parametrize physics beyond the Standard Model. 
The optimistic interpretation of all this is that, just like Fermi anticipated particle physics at 
100 GeV in 1934, the electroweak gauge theory requires new physics to tame the divergences 
associated with the Higgs potential. By the most conservative estimates this new physics 
is within our reach. Avoiding fine tuning requires A < 2 — 3 TeV to be revealed by the 
CERN Large Hadron Collider (LHC), possibly by the Tevatron. For instance, for M# = 
115- 200 GeV, 

Sv 2 

— <10^A = 2- 3 TeV . (32) 
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Any physics which may turn on beyond the electroweak scale can easily enhance weak 
interaction cross sections beyond Standard Model predictions. In this context, neutrino 
telescopes provide an important probe of new ideas in particle physics. One caveat for the 
realization of such a program is that the event rates depend on the as yet unknown extrater- 
restrial neutrino flux, generating an ambiguity in obtaining cross sections from event rates. 
However, the event rates for up- coming and down-going neutrinos have different responses 
to the inelastic cross section [6jJ. Therefore, a full study in the two-dimensional parameter 
space (cross section and flux) can serve as a powerful testbed for probing anomalous neu- 
trino interactions |62l . [63[ . Bounds (or discovery potential) of neutrino fluxes emerge in a 
correlated fashion with the cross section probes. 

Physics at the high energy frontier is the physics of partons. We master this physics 
with unforeseen precision because of a decade of steadily improving HERA measurements 
of the nucleon structure These now includes experiments using targets of polarized 

protons and neutrons. HERA is our nucleon microscope, tunable by the wavelength and 
the fluctuation time of the virtual photon exchanged in the electron proton collision. The 
wavelength of the virtual photons probing the nucleon is reduced with increased momentum 
transfer Q. The proton has now been probed to 10~ 3 fm, about one-thousandth of its size. 
In the interaction, the fluctuations of the virtual photons survive over distances ct ~ 1/x, 
where x is the relative momentum of the parton. HERA now studies the production of chains 
of gluons as long as 10 fm, an order of magnitude larger than and probably insensitive 
to the proton target. These are novel QCD structures, the understanding of which has 
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been challenging. The x — Q 2 region probed by cosmic neutrinos lies in this challenging 
region, beyond the reach of the HERA experiments. It is therefore neccesary to think of 
reliable methods for extrapolating the structure functions to very low x. For the purpose of 
phenomenological studies, we have chosen to use the parametrizations of the CC 

/ p \ 0.363 

acc = 5 - 53 lGeV) Pb ' (33) 



0.363 



and NC 



neutrino-nucleon cross sections, as derived in |65( using the CTEQ4 parton distribution 
functions |6^ |. 

Deviations of the neutrino-nucleon cross sections from the simple parton model can result 
from more sophisticated treatments of QCD or from new physics beyond the Standard Model. 
Included in the former are saturation effects which can substantially modify simple parton 
results at small x [67] . These effects can significantly alter the total cross section at high 
energies |68j], softening the power law behavior predicted by the simple parton model toward 
compliance with the Froissart bound . At the forefront of the new physics category are 
scenarios based on TeV-scale quantum gravity [7(|. Enhancements of the neutrino cross 
section in this framework may come from exchange of towers of Kaluza- Klein gravitons [z]J 
and black hole production [72]. Similar enhancements may originate from TeV-scale string 
excitations [73|, electroweak instanton processes jlll, or in some supersymmetric models 
through direct channel production of superpartner resonances [T5| . 

In what follows we discuss the potential of IceCube to reveal deviations of the neutrino- 
nucleon cross section from Standard Model predictions or to set limits on such deviations, 
without assuming particular neutrino fluxes or cross sections. It will be convenient, however, 
to present results relative to standard reference values. For the reference cross section, we 
adopt &sm the neutrino-nucleon cross section of the simple parton model given in Eq. (|3~3*)). 2 
For the reference flux, we adopt the WB flux given in Eq. (|25[) . Altogether, in the energy 
range 10 7 ~ 75 GeV, after appling selection criteria neccesary for good energy resolution, we 
estimate an effective integrated luminosity (over the lifetime of the experiment) £ ~ 30 nb -1 . 
At an average energy (E) = 10 7 ' 25 GeV, ~ 4 x 1CT 23 GeV -1 cm -2 s _1 sr _1 per flavor 
and (Tsm({E)) ~2x 1CT 33 cm 2 . In order to separate the flux and cross section contributions 
to the event rate, we require two independent pieces of data, which hereafter we identify as 
up-coming and down-going events. 

For down-going events, the probability of neutrino conversion is always small, barring 
extraordinary enhancements to neutrino cross sections. Letting o v n denote the total (NC 
+ CC) neutrino-nucleon cross section, the down-going event rate is therefore 

\f - r ^ ° vN 

■'vdown — '-'down , v > yJO ) 

<PWB a SM 

where <p u indicates the average extraterrestrial v + V flux per flavor in the energy bin of 
interest. The constant Cd OW n depends on exposure and varies according to neutrino flavor 



2 We neglect NC interactions, which at these energies serve only to reduce the neutrino energy by approxi- 
mately 20%. 
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from experiment to experiment. Here, a V N indicates any enhancement of the cross section 
which will increase the event rate for down-going neutrinos, but because of absorption, also 
suppress the up-coming events. The latter can be achieved through cuts on the concommi- 
tant increase in shower energy fraction greater than or equal to that characterizing the CC 
Standard Model process. 

For up-going events, the dependence on cross section is completely different. At energies 
above 10 6 GeV, neutrino interaction lengths become smaller than the radius of the Earth, 
Rq ~ 6371 km. This implies that most upward-going neutrinos are shadowed by the Earth, 
and only those that are traveling at large nadir angles along chords, t = 2R S) cos8, with 
lengths of order their interaction length, l v (E u ) = [p CT ust-^4 &v n{ E v ) \~ 1 , can produce a visible 
signal (p cr ust — 2.65 gem -3 is the density of the Earth's crust) [76j. 

The dependence of upward-going event rates on anomalous neutrino cross sections de- 
pends on the source of the anomaly. We consider two prominent cases. First, in many new 
physics cases, the Standard Model CC cross section <7sm is not altered, but there are new 
neutrino interactions that produce showers. For E v > 10 7 GeV, the neutrino interaction 
length satisfies li -C l v < R®, and so the up-going event rate is [62| 

^ P = C up ^ldM (^>l) , (36) 



where l\ is the typical lepton path length in Earth. Extreme enhancements to <j v n may 
reduce l v to //, leading to a different parametric dependence in Eq. (j3BJ), but such cases are 
neglected here. Next we consider the possibility of screening. In order to probe deviations 
from the unscreened parton model calculations of the cross section, it is necessary to note 
that the screening correction affects CC and NC equally. In this case one finds |6l|, |62j , 

\f _ -^screen T* ^SM fO r 7\ 
•A/ U p - O up — — , [6() 

where <t u n and o"sm are CC cross sections with and without screening, respectively. 

Constraints on the two-parameter space can be obtained using AMANDA data. To see 
this, we note that the 90% CL upper bound on the diffuse neutrino flux 

E 2 V <f) u max = 3.3 x 1(T 7 GeV cm" 2 s" 1 sr" 1 (38) 



(per flavor) reported by the AMANDA Collaboration |36J was derived on the basis of CC 
neutrino interactions within Standard Model physics, assuming an E~ 2 dependence of the 
flux, valid for 10 6 GeV to 10 9 5 GeV neutrinos. In full generality, the stated quantity 
max constitutes two joint constraints on two unknowns: the diffuse neutrino flux and any 
deviation of the simple parton model neutrino-nucleon cross sections. Since the energy 
distribution of the AMANDA data peaks in the energy bin of interest, it is reasonable to 
use max ((-E)) as the upper limit in the bin. Applied to down-going events, the constraint 
implies 

" " < Cax • (39) 



Dividing Eq. (j3~9j) by 4>wb gi ves 



C"SM 

^r- — < 26 (40) 
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at 90% CL. A similar analysis for up-going events yields 




< 26 



(41) 



for the case of new physics contributions, and 



<P U o" SM 



< 26 



(42) 



for the case of screening. These constraints exclude the shaded and cross-hatched regions of 
Fig. [7| The shaded region is excluded by down-going data, and the cross-hatched region is 
excluded by the up-going data, assuming screening. The shaded and cross-hatched regions 
meet at ovjv/osm = 1- As a result, irrespective of cross section assumptions, one finds an 
upper bound on the neutrino flux in the energy range 10 7 GeV to 10 7,5 GeV of <f) u < 26(p^ VB . 

How will these results improve in the near future? To evaluate the prospects for IceCube, 
one must determine Cd OW n and C up . We focus on the case in which new neutrino interactions 
modify the NC cross sections, but leave the CC cross sections invariant. In NC processes 
most of the energy is carried off by pions. At TeV energies, the interaction mean free path of 
tt^ in ice is orders of magnitude smaller than the pion decay length, and so nearly all energy 
is channeled into electromagnetic modes through ir° decay. To estimate the efficiencies for 
down-going events in the NC channel, we therefore adopt as our basis of comparison the 
electromagnetic showers induced by v e . With this in mind, 



where tit((E)) ~ 6 x 10 38 . Inserting these numbers into Eq. ()43|) . along with a lifetime 
of the experiment T = 15 yr and AE = 2.2 x 10 7 GeV, one obtains Cd OW n — 3. The 
v e atmospheric background is taken to be zero, (the potential "prompt" component can be 
eliminated by taking advantage of the isotropic ^-channel) . To determine the corresponding 
quantity for up-going events C up , we first note that the absence of oscillation precludes a 
v T atmospheric background. For this reason the detection prospects for up-going neutrinos 
(with 10 7 GeV < E v < 10 7 5 GeV) are brighter for v T than the other flavors, and we focus 
on them below. From Figs. Q and 01 after correcting for the Earth's absorption effects, one 
obtains the v T effective aperture, (A£l) c g ~ n km 2 sr. Thus, the normalization constant for 
r's showering in IceCube is found to be 



down 



27r T J 0wb a SM n T dE u 
2ir<% m ((E)) a SM ((E)) n T ((E)) TAE , 



(43) 




(44) 



where [Z2| 




J 2x 1(T 3 0^ B (£„) dE, 
5.5 x 1CT 1 km~ 2 yr -1 sr~ 



V 



-1 



(45) 



is the r-lepton flux produced in v T interactions inside the Earth for a v ^ = osm- 
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FIG. 7: Projected determinations of neutrino fluxes and cross sections at ~ 6 TeV from future 
IceCube data. Two cases with (A/J^ n , A/"°p S ) = (3,25) and (50,25) are considered. In both cases, 
the best fit flux and cross sections are shown, along with the 90%, 99%, and 99.9% CL exclusion 
contours. Contours of constant A/" up = 25 (dotted), A/down = 3 (left dashed), and A/down = 50 (right 
dashed) are also shown. Neutrino fluxes and cross sections excluded by AMANDA at 90%CL [3^| 
are also indicated: The shaded (cross-hatched) region is excluded by null results for down-going 
(up-going) events psT ]. 



Given these estimates of Cdown and C np , one can now determine projected sensitivities of 
IceCube to neutrino fluxes and cross sections. The quantities Af up and A/down as defined in 
Eqs. and (pUJj) can be regarded as the theoretical values of these events, corresponding 
to different points in the a u n / &sm, 0^/0wb parameter space. For a given set of observed 
rates A/"°p S and A/J^, two curves are obtained in the two-dimensional parameter space by 
setting A/°p S = A/up and A/Jown = A/down- These curves intersect at a point, yielding the most 
probable values of flux and cross section for the given observations. Fluctuations about 
this point define contours of constant \ 2 m an approximation to a multi-Poisson likelihood 
analysis. The contours are defined by 

x 2 = J2 2 Wi ~ ^ ohs ] + 2 K pbs ln [K obs /^i] , (46) 

i 

where i = up, down [78]. 

In Fig. we show results for two representative cases that are consistent with the 
AMANDA bounds derived above. In the first case, we assume a u jy = <tsm and (p u = 0wb> 
leading to 3 down-going and 25 up-going events. The 90%, 99%, and 99.9% CL contours 
are those given in the lower left of the figure. (These contours will be slightly distorted for 
Vvn/vsm < 1, where Eq. (}3l)j) receives corrections, but we neglect this effect.) We see that, 
even in the case that event rates are in accord with standard assumptions, the neutrino- 
nucleon cross section is bounded to be within 40% of the Standard Model prediction at 90% 
CL. This is at a center-of-mass energy y/s ~ 6 TeV, far beyond the reach of any future 
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man-made accelerator. 3 In the second case, we consider a scenario in which the number of 
observed upcoming events remains at 25, but the number of down-going events is 50. In the 
second case, clearly one has discovered new physics at well beyond ha. 

The analysis discussed in this section is very general as the technique to extract the 
cross section is model independent and can be readily extended to probe higher energy bins. 
From Fig. EJone can see that, if there are hidden sources in the heavens, we expect Cd OW n 
and C up to be about a factor of 15 higher, making it possible for IceCube to probe 40% 
(70%) enhancements from Standard Model predictions at the 90% (99.9%) CL after only 1 
year of observation. 



V. FLAVOR METAMORPHOSIS 
A. Neutrino Oscillations 

The Standard Model is based on the gauge group 

G?M ge = SU(3) C x SU(2) L x U(1) Y , (47) 

with three fermion generations. A single generation consists of five different representations 
of the gauge group: Q L (3, 2, 1/6), U R (3, 1, 2/3), D R (3, 1, -1/3), L L (1, 2, -1/2), E R (1 1, -1); 
where the numbers in parenthesis represent the corresponding charges under G|^ ge . The 
left handed fermions belong to SUl{2) doublets while the right handed chiralities to SUl{2) 
singlets. The model contains a single Higgs boson doublet, #(1,2,1/2), whose vacuum 
expectation value breaks the gauge symmetry Gsm SG SU(3)c x [7(1)em- The Standard 
Model also comprises an accidental global symmetry 

Gf° bal = U(1) B x U(l) e x U(l)p x U(1) T , (48) 

where U(1)b is the baryon number symmetry, and £/(l) ej/v - are three lepton flavor symme- 
tries, with total lepton number given by L = L e + L M + L T . It is an accidental symmetry 
because we do not impose it. It is a consequence of the gauge symmetries and the low en- 
ergy particle content. It is possible (but not necessary), however, that effective interaction 
operators induced by the high energy content of the underlying theory may violate sectors 
of the global symmetry. 

In the Standard Model charged lepton (and also quark) masses arise from Yukawa inter- 
actions, which couple the right-handed fermion singlets to the left-handed fermion doublets 
and the Higgs doublet 

£^ awa = \iL L H E R + h.c. (49) 

After spontaneous electroweak symmetry breaking these interactions lead to charged fermion 
masses, mi = \i v/y/2, but leave the neutrinos massless. No Yukawa interaction can be 
written that would give a tree level mass to the neutrino because no right-handed neutrino 



3 For small x physics, Standard Model predictions for neutrino-nucleon collisions at ^J~s ~ 6 TeV probe x 
down to Myy/s ~ 2 x 10~ 4 79] ; LHC will probe parton center-of-mass energies \/S > 50 GeV (twice 
the minimum jet trigger threshold), corresponding to x > 50/14000 ~3x 10~ 3 . A future hadron collider 
would require y/s ~ 200 TeV to sample ik2x 10 -4 . 
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field exists in the Standard Model. One might think that neutrino masses could arise from 
loop corrections. This, however, cannot be the case, because any neutrino mass term that can 
be constructed with Standard Model fields alone would violate the global lepton symmetry. 
In other words, the Standard Model predicts that neutrinos are strictly massless, and thus 
there is neither mixing nor CP violation in the leptonic sector. 

In recent years, stronger and stronger experimental evidence for nonzero neutrino masses 
has been accumulating 80]. The weak interaction coupling the W boson to a charged lepton 
and a neutrino, can also couple any charged lepton mass eigenstate l a (a = e, p, r) to any 
neutrino mass eigenstate Uj (j = 1,2,3) . The superposition of neutrino mass eigenstates 
produced in association with the charged lepton of flavor a, \u a ) = U*j\uj), is the state 
we refer to as the neutrino of flavor a, where U a j's are elements of the neutrino mass-to- 
flavor mixing matrix, called the Maki-Nagakawa-Sakita (MNS) matrix. The density matrix 
of a flavor state, p a , can be expressed in terms of mass eigenstates by p a = \v a )(v a \ = 
YlijUaiUajWilivjl- This is a pure quantum system, therefore the density matrix satisfies 
Tr p 2 = Tr p = 1. The time evolution of the density matrix 

% = -i [H, P] , (50) 

is governed by the Hamiltonian H of the system, thereby after traveling a distance L an 
initial state u a becomes a superposition of all flavors, with probability of transition to flavor 
13 given by P Uot ^ = Tr\p a (t)p^\, or equivalently 

Pv^vp = <W - 4 ^ e (Ki Upi U aj U* pj ) sin 2 A 4j 

i>j 

+ 2 Sm (Ki Upi U aj U; 3 ) sin 2Ay , (51) 

i>j 

where A^- = 1.27 Sm^L/ E v , when L is measured in km, E v in GeV and Sm 2 j = m 2 — m 2 in 
eV 2 . 

The simplest and most direct interpretation of the atmospheric data [81] is that of muon 
neutrino oscillations. The evidence of atmospheric disappearing is now at > 15<r, most 
likely converting to v T . The angular distribution of contained events shows that for E v ~ 
1 GeV, the deficit comes mainly from L atm ~ 10 2 — 10 4 km. The corresponding oscillation 
phase must be maximal, A atm ~ 1, which requires Sml tm ~ 10" 4 - 10~ 2 eV 2 . Moreover, 
assuming that all upgoing z/^'s which would yield multi-GeV events oscillate into a different 
flavor while none of the downgoing ones do, the observed up-down asymmetry leads to a 
mixing angle very close to maximal, sin 2 26* atm > 0.85. These results have been confirmed 
by the KEK-to-Kamioka (K2K) experiment which observes the disappearance of accelerator 
v^s at a distance of 250 km and find a distortion of their energy spectrum with a CL of 
2.5 — 4cj The combined analysis leads to a best fit-point and la ranges, Sm 2 tm = 

2.2t° i x 10" 3 eV 2 and tan 2 # atm = lt° H 0- On the other hand, reactor data requires 
| Z7e3 1 2 1. Thus, to simplify the discussion hereafter we use the fact that \U e 3\ 2 is nearly zero 
to ignore possible CP violation and assume that the elements of the MNS matrix are real. 
The simple vacuum relations between the neutrino states produced by the weak interactions 
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are then 



1 
1 



ki> = -^sinfl©^) - \v T )) - cos6 Q \v e ) 



\v 2 ) = -j= cos Oq^) - \u T )) + sin 9 Q \u e ) (52) 



1 

3" 



1^3) = + Wr)) 



where 8 p i s the solar mixing angle. Data collected by the Sudbur y N eutrino Observatory 
(SNO) [84j in conjuction with data from SuperKamiokande (SK) |4l| show that solar v' e s 
convert to or v T with CL of more than la. The KamLAND Collaboration [85] has 
measured the flux of V e from distant reactors and find that TVs disappear over distances 
of about 180 km. The combined analysis of Solar neutrino data and KamLAND data are 
consistent at the 3<r CL, with best-fit point and la ranges: 5m 2 & = 8.2^3 x 1CT 5 eV 2 and 
tan 2 # = 0.3918$ B- 

Altogether, neutrinos are massive and therefore the Standard Model needs to be extended. 
The most economic way to get massive neutrinos would be to introduce the right-handed 
neutrino states (having no gauge interactions, these sterile states would be essentially un- 
detectable and unproduceable) and obtain a Dirac mass term through a Yukawa coupling. 
However, in order to get reasonable neutrino masses (below the eV range) the Yukawa cou- 
pling would have to be unnaturally small \ v < 10~ n (note that for charged fermions the 
Yukawa couplings range from X t ~ 1 for the top quark down to A e ~ 10~ 5 for the electron). 
Other attempts to give masses to neutrinos require extension of the Standard Model in 
more radical ways. Throughout this review we will not concern ourselves with the origin 
of the neutrino masses (for a comprehensive discussion of this subject see e.g., j!3|), but 
simply assume that neutrinos acquire mass by some mechanism and discuss the reach of 
IceCube in probing physics beyond conventional neutrino mass-induced oscillations. The 
methods discussed in this section are readily applicable to speculations concerning other 
non-conventional physics associated with neutrino flavor mixing. 

The only interaction a neutrino can feel, besides the weak one, is the gravitational interac- 
tion. A gravitational field, could contribute to neutrino oscillations if the different neutrino 
flavors are affected differently by gravity, i.e., if the Equivalence Principle is violated |88j. 
The deviation can be parametrized by assuming that the post-Newtonian parameters in the 
expansion of a given metric have a species dependent value, representing thus the effective 
coupling constant of the different flavors to a given geometry. In the weak field limit of a 
static source, to first order in the Newtonian potential <ft — GM/r <C 1, one can consider 
the metric to be g tt — 1 — 2-f(fi and = — %(1 + 270) . The parameters 7 and 7 are theory 
dependent (e.g., 7 = 7 = 1 in General Relativity), but for any given theory, their value is 
universal (the same for all kinds of particles) only if the Principle of Equivalence is satisfied. 
Otherwise their value could be different for test particles with different internal quantum 
numbers and or different energies. 

Of particular interest here, it has been shown [89( that the value of 7 is the same to an 
accuracy of about 0.1% for neutrinos and photons received from the SN87a. Moreover, if 
neutrinos are light, 7 is the same for neutrinos (and antineutrinos |9(|) of different energies 
(ranging from 7 to 40 MeV) up to an accuracy of one part in 10 6 . It is worthwhile to note 
that the above results have been obtained with the hypothesis 7 = 1. This hypothesis holds 
with a data set containing only one type of particle, as possible deviations of 7 from 1 can be 
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absorbed by redefining the mass of the source. The hypothesis may not be justified however, 
if we compare different particles or particles with different energies, since the value of 7 
can be particle dependent or energy dependent just like 7. The laboratory experiments on 
neutrino oscillations testing the universality of the gravitational redshift provide information 
on the possibility that the value of 7 be different for different neutrino flavors. If there are 
deviation from universality in the value of 7, in fact, the energy splitting induced by gravity 
may contribute in general to the transition probability between different flavors, even if 
gravity is not the primary source of oscillations, i.e., even if the energy eigenstates are 
unchanged by the weak gravitational field present in the laboratory. 

At this stage it is worthwhile to point out that for a constant potential 0, a violation 
of the equivalence principle is phenomenologically equivalent to the breakdown of Lorentz 
invariance resulting from different asymptotic values of the velocity of the neutrinos, c\ 7^ C2, 
with vx and v 2 being related to and v T by a rotation of an angle £ v ij [9l[ . 

Gravity induced neutrino oscillations can also arise if the spacetime manifold has a non- 
symmetric Christoffel connection j^. This is because the flavor eigenstates produce in the 
charged weak interactions are not, in general, the eigenstates of the torsion Hamiltonian, 
and then the oscillations can be strongly damped or enhanced in the presence of strong 
gravitational fields. The torsion eigenstates can be interpreted as different mixtures of 
right-handed and left-handed fields, yielding a deviation of the weak neutrino current from 
the standard V — A form. For this vector-like interaction the gravity-induced oscillation 
wavelength is energy independent. Violation of CPT resulting from Lorentz- violating effects 
also leads to an energy independent contribution to the oscillation wavelength [93| which is 
a function of the eigenvalues of the Lorentz violating CPT-odd operator and the rotation 
angle between the corresponding neutrino eigenstates z/j and the flavor eigenstates v a . 

In the most general situation in which neutrinos are massive and gravity may contribute 
to the transition probability between different flavors, we may have \v a ) 7^ \vi) 7^ \pg)i 
where \vq) denotes hypothetical gravity eigenstates. The eigenstates of the total energy 
are obtained by diagonalizing the matrix which includes the contributions of the mass and 
gravitational terms. In the case of 2 flavor mixing, the total matrix for neutrino (+) or 
antineutrino (— ) oscillations can be written as 

where accounts for a possible relative sign of the gravity effects between neutrinos and 
antineutrinos 4 and S n parametrizes the size of the new physics terms. The matrices and 
U €n ,±^ are given by 

TT _ ( cos 9 sin6»\ / cos£ n sin£ n e ±M7n \ , 

Ue \-sin6 cos 6 J ' U €».±»J» ^- S in^e^ cos^n J ' 1 ' 

where i] n denotes the possible non-vanishing relative phases. Note that in contrast to the 
conventional oscillation length, new physics predicts neutrino oscillations with wavelengths 
that are either constant or decrease with energy. Therefore, IceCube with an energy reach 
in the 10 2 GeV < E v < 10 6 GeV range for atmospheric neutrinos, is the ideal experiment 



4 For the CPT violating case <r + = —a , otherwise a + = a 
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to search for new physics. Moreover, these neutrinos have high enough energy for the 
standard oscillations (due to Sm^) to be very suppressed so the z/s behave essentially as 
flavor eigenstates. For simplicity in what follows we concentrate on oscillations resulting 
from tensor-like interactions that lead to an oscillation wavelength inversely proportional to 
the neutrino energy. The results can be directly applied to oscillations due to violations of 
the equivalence principle with the identification, £ x = £ vep and 5i = 2|0|(7x — 72) = 2|0| A-y, 
as well as oscillations due to violations of local Lorentz invariance with the identification, 
£1 = £vii and 5i = Ac/ c. 

If the effects of new physics are constant along the neutrino trajectory, the survival 
probability for total oscillation length L to t takes the form 



Pu^u, = 1 - Pv^ut = 1 - sin 2 20 sin 2 (Ay K) , (55) 
where 1Z = L atm /Ltot; ^ the contributions of mass and gravity are both non- vanishing then 
n = \/l + R 2 + 2R (cos 2# atm cos 2f 1 + sin 26 atm sin 2^ cos rj{j . (56) 
Here, B is the rotation angle which diagonalizes the matrix in Eq. (|53j) . 

sin 20 = — - (sin 26 atm + R 2 sin 2 2f a + 2R sin 26 atm sin 2f 1 cos?7i) , {01) 
1Z 2 



and 



T-, -^atm 



where Lq defines the gravity oscillation length. 

The experimental probes of neutrino oscillations, measuring P, provide simultaneously 
information on and L tot . In the absence of positive results, however, an upper bound on 
P does not fix any value for G and L tot separately: it determines only an allow region in 
the (0,L tot ) plane, or equivalently in the plane. At present, the strongest limits on 

neutrino oscillations due to new physics come from the non-observation of departures from 
the standard oscillations seen in atmospheric neutrinos at SK, and the confirmation of 
oscillations with the same oscillation parameters (i.e., ^^ 2 tm) from K2K. The 3a bounds 
from the updated combined analysis [95j of SK and K2K data are shown in Fig. |HJ 

For most of the neutrino energies considered here, the standard mass-induced oscillations 
are suppressed and the gravity effect is directly observed. As a consequence, the results will 
be independent of the phase 7/1 and one can set the new physics parameters to be in the 
range 5\ > and < £1 < tc/4. 

The Hamiltonian of Eq. ()53|) describes the coherent evolution of the v li —v T ensemble 
for any neutrino energy. High-energy neutrinos propagating in the Earth can also interact 
inelastically with the Earth matter and as a consequence the neutrino flux is attenuated. For 
atmospheric neutrinos oscillation, attenuation, and regeneration effects occur simultaneously 
when the neutrino beam travels across the Earth's matter. For the phenomenological analysis 
of conventional neutrino oscillations this fact can be ignored because the neutrino energies 
covered by current experiments are low enough for attenuation and regeneration effects to be 
negligible. Especially for non-standard scenario oscillations, future experiments probe high- 
energy neutrinos for which the attenuation and regeneration effects have to be accounted 
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FIG. 8: Sensitivity limits in the (Ac/c,£ v ii) and (2\<p\ A7, £ vep ) planes at 90%, 95%, 99% and 
3<T CL. To account for the uncertainty associated with the poorly known prompt neutrino fluxes 
results for two extreme models are shown: TIG [la ] (full lines) and RQPM [2l( (filled regions). 
The difference is about 50% in the strongest bound on 5\. The hatched area in the upper right 
corner is the present 3a bound from the analysis of SK and K2K data 0]. 



for simultaneously. Attenuation effects due to CC and NC interactions can be introduced 
by relaxing the condition Tr(p) = 1, so that Eq (joTJj) can be re- written as 

!£ = -i[H,p]- y £±.{n a ,p}, (59) 

a u 

where l u is the neutrino mean free path and Tl a = v a ® i/t is the u a state projector. 

Neutrino oscillations introduced by new physics effects result in an energy dependent 
distortion of the zenith angle distribution of atmospheric muon events. This effect can be 
quantify by evaluating the expected angular and E^ n distributions in the IceCube detector 
using Eq. © in conjunction with (and v^) fluxes obtained after evolution in the Earth for 
different sets of new physics oscillation parameters. Together with i^-mduced muon events, 
oscillations also generate \x events from the CC interactions of the v T flux which reaches 
the detector producing a r that subsequently decays as r — > iw^v T and produces a fi in 
the detector. Generalizing the formulae given in Sec. |TT] it is straightforward to see that the 
number of i^-induced muon events is 

/I poo roc poo poo poo 

dcosO / dl^ / dEl n / dEl / dE T / dE v (60) 
1 Jl™ in Jm^ J JE° J E T 

dFVT f TTi „ _ _ m rfcr CC , rn m \ . - dN dcc „ Q fin Q 



-(E u , cos6)^±(E u , E T ) n T ^(E T , E°)E(E° K fin , /„) A c c 



dEu^e" v ' j ^eV" v ' TJ 1 dEj ^' ^ '^~ cff ' 
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RPQM 
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log 10 (^ n /GeV) 


-1 < cos6> < 


-0.6 -0.6 < cos# < -0.2 


-1 < cos# < 
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2.00 
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24647 
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3.20 
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1289 
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- 4.00 
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410 
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498 
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96 
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- 4.60 


86 
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36 


89 
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- 4.80 


38 


74 


13 


36 


4.80 


- 5.00 


16 


34 


5 


14 


5.00 


- 9.00 


10 


28 
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TABLE I: Number of expected atmospheric z^-induced muon events in 10 years of IceCube oper- 
aton in the 16 energy bins and angular bins discuss in the text, assuming no new physics effect is 
observed. 



where dN dec /dE® is given in Ref. [96(. To quantify the energy- dependent angular distortion 
one defines the vertical-to-horizontal double ratio 



R 



hor / ver 



fin.i 



P, 



hor 



fin,-. 



Af 5 ^\El n '\ -0.6 < cosfl < -0.2) 
A/* 1=0 (£j n '\ -0.6 < cosfl < -0.2) 
-1< costf < -0.6) 



(61) 



\Ep\-l < cos6 < -0.6) 

where -E* n ' 4 denotes integration in an energy bin of width 0.2 log 10 (^ n ' i ). This leads to 
16 Ej^ bins in the data sample: 15 bins between 10 2 and 10 5 GeV and one containing all 
events above 10 5 GeV. Note that the double ratio eliminates uncertainties associated with 
the overall normalization of the atmospheric fluxes at high energies. Additionally, in the 
definition of the double ratio we have conservatively included only events well below the hori- 
zon cos 9 < —0.2 to avoid the possible contamination from missreconstructed atmospheric 
muons which can still survive after level 2 cuts in the angular bins closer to the horizon. In 
order to estimate the expected sensitivity we assume that no new physics effect is observed 
and define a simple X 2 (^i)£i) function as 



16 

E 

i=l 



(Rhor/ v er{Ef i n,t ) 



1) 



cr, 



(62) 



stat.j 



where cr sta t,i is computed from the expected number of events in the absence of new physics 
effects. The expected event rates are given in Table |TJ 
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Figure |H1 shows the sensitivity limits in the [Si, £i]-plane at 90%, 95%, 99% and 3cr CL 
obtained from the condition x 2 (3ii£i) < Xmax(CL, 2dof). The figure illustrates the improve- 
ment on the present bounds by more than two orders of magnitude even within the context 
of this very conservative analysis. The loss of sensitivity at large Si is a consequence of the 
use of a double ratio as an observable. Such an observable is insensitive to new physics 
effects if Si is large enough for the oscillations to be always averaged leading only to an 
overall suppression. 



B. Quantum Decoherence 

What is the meaning of quantum gravity? It means that spacetime itself is subject to 
quantum laws, necessitating inherent fluctuations in the fabric (metric and topology) of space 
and time. These microscopic boiling bubbles force on spacetime a foam-like structure (9?}. 
A heuristic example pictures spacetime to be filled with tiny virtual black holes that pop in 
and out of existence on a timescale allowed by Heisenberg's uncertainty principle [98(. These 
black holes conserve energy, angular momentum, and electric and color charge (unbroken 
gauged quantum numbers), but they are believed not to conserve global quantum numbers. 
The transition between initial and final density matrices associated with black hole formation 
and evaporation is not factorizable into products of ^-matrix elements and their hermitian 
conjugates. The evolution of such a quantum system is characterized by a superscattering 
operator § that maps initial mixed states to final mixed states, p out = §Pm, with S ^ 
S^S [99]. In other words, there may be a loss of quantum information across the black hole 
event horizons, provi ding an environment that can induce decoherence of apparently isolated 
matter systems |l00j ]. 5 

Of particular interest is the question whe ther fast baryon decay can proceed via virtual 



black hole states in the spacetime foam |101| . The process is envisioned as the simultaneous 
absorption of two quarks into the black hole, followed without memory of the initial state by 
the thermal emission of an antiquark and a lepton, entailing a change in the global baryon 
and lepton quantum numbers qq — > ql. The probability that two quarks in a proton of size 
Aq^ d pass within a fundamental Planck length, within the Heisenberg lifetime uncertainty 
of the black hole is oc (Aqcd/M^g) 4 , where Mqq is th e gravitational UV cutoff. Thus , the 
present limit on the proton lifetime, t p ~ 10 33 yr 102], implies Mqq > 10 16 GeV |l0ll |. 6 

Measurements of flavor transformations in a neutrino beam can provide a clean and 
sensitive probe of interactions with the spacetime foam. Without interference from the 
gravitational sector, oscillations in the neutrino sector provide a pure quantum phenomenon 
in which the density matrix has the properties of a projection operator. Because black holes 
do conserve energy, angular momentum (helicity), color and electric charge, any neutrino 
interacting with the virtual black holes needs to re-emerge as a neutrino. As an example, 
if spacetime foam black holes do not conserve U(l) e x U(l)^ x U(1) T , neutrino flavor is 



5 Recent developments in string theory show that there is in fact no information loss in the fundamental 
underlying theory, nevertheless such information loss might appear in the effective low energy theory. 
This is clearly an important issue, worthy of investigation. 

6 Note that interactions through the higher dimensional QQQL operator can be prevented if one separates 
the quark and lepton fie lds f ar enough in an extra dimension, so that their wave function overlap is 
exponentially suppressed [l03j . 
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randomized by interactions with the virtual black holes. The result of many interactions is 
to equally populate all three possible flavors. 

The Hamiltonian evolution in Eq. (|5U|) is a characteristic of physical systems isolated from 
their surroundings. The time evolution of such a quantum system is given by the continuous 
group of unitarity transformations, U t = e~ lHt , where t is the time. The existence of the 
inverse of the infinitesimal generator, H, which is a consequence of the algebraic structure of 
the group, guarantees reversibility of the processes. For open quantum-mechanical systems, 
the introduction of dissipative effects lead to modifications of Eq. (joT?|) that account for the 
irreversible nature of the evolution. The transformations responsible for the time evolution 
of these sys tems are defined by the operators of the Lindblad quantum dynamical semi- 
groups [104]. Since this does not admit an inverse, such a family of transformations has the 
property of being only forward in time. 

The Lindblad approach to decoherence does not require any detailed knowledge of the 
environment. The corresponding time evolution equation for the density matrix takes the 
form 

^ = - i [H cS ,p}+V[p} , (63) 
where the decoherence term is given by 

^] = 4E([ & ^ & !] + [^' & !]) • ( 64 ) 

i 

Here, H e ^ = H + H^ is the effective Hamiltonian of the system, H is its free Hamiltonian, 
accounts for possible additional dissipative contributions that can be put in the Hamiltonian 
form, and {bj} is a sequence of bounded operators acting on the Hilbert space of the open 
quantum system, 7i, and satisfying Yljbjbj G B(H), where B(H) indicates the space of 
bounded operators acting on H. 

The intrinsic coupling of a microscopic system to the space-time foam can then be in- 
terpreted as the existence of an arrow of time which in turn makes possible the connection 
with thermodynamics via an entropy. The monotonic increase of the von Neumann entr opy 
S(p) = — Tr (p hip), implies the hermiticity of the Lindblad operators, bj = b^ jl05j | . In 
addition, t he c onservation of the average value of the energy can be enforced by taking 
JLbj] -0 HH}. 

In a 2-level system with mixing angle 9 (relevant f or a tmospheric neutrinos), the 
survival probability after propagation of a distance L is 

Pu^u, = 1 - Pu^u T = 1 - sin 2 29 sin^e-^Atf] , (65) 

where 7 has dimension of energy, and its inverse defines the typical (coherence) length after 
which the system gets mixed. 7 Thus, for 7L ~ 0(1) one expects significant deviations from 
standard mass- induced oscillations formula P u „-«/„ = 1 — sin 2 29 sin 2 A^. For dissipative 



7 The parametrization that lead s to Eq. I|65[) is appropriate for including dis sipat ion as a perturbation on 
the standard oscillations |l0sj . In the most general solution of Eq. (|63(l |l09| . after propagation over 
extreme long-baselines (i.e., L — > 00) the system evolves to an equal flavor ratio regardless of the initial 
flavor content and mixing angle. 
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scenarios in which decoherence effects vanish in the weak gravitational limit, Mqg — > oo, 7 
can be parametrized as 

^G^)"®""^' (66) 

where k is a dimensionless parameter, which by naturalness is expected to be 0(1), and 
n > 2. 

The energy behavior of 7 depends on the dimensionality of the operators bj. But care 
must be taken, since T> is bilinear in the bj, and due to the hermiticity requirement, each bj 
is itself at least bilinear in the neutrino fields ip. Examples are 

bj oc / d 3 x i) ] {id t ) j ip , (67) 



where j = 0, 1, 2 .... A Fourier expansion of the fields ip, inserted into Eq. ()67|1. gives 
the energy behavior bj cx E 3 V , and hence 7 oc E%>. This restriction of the energy behavior 
to non-negative even powers of E u may possibly be relaxed when the dissipative term is 
directly calculated in the most general space-time foam background. 

An interesting example is the case where the dissipative term is dominated by the 
dimension-4 operator bi, J d 3 x ip' id t i/j , yielding the energy dependence 7 cx E^/Mqq. This 
is characteristic of non-critical string theories where the space-time defects of the quan- 
tum gravitational "environment" are take n as recoiling Z)-branes, which generate a cellular 
structure in the space-time manifold lid! . 



A best fit to data collected by the SK atmospheric neutrino experiment |41|, allowing for 
both oscillation and decoherence yields, for n = 2, 

* (^)" <0 - 9Xl °" 7 ' (68) 

at the 90% CL |l07j |. A similar lower limit on Mqg has been obtained at Fermilab. The 
CCFR detector is sensitive to <-> u e 111 ] and u e <->• v T 112l | flavor transitions. Neutrino 



energies range from 30 to 600 GeV with a mean of 140 GeV, and their flight lengths vary 
from 0.9 to 1.4 km. A best fit to the data allowing for both oscillation and decoherence 
yields 

' M QgV X .on., in-24 



GeV 



< 2.0 x 10" i4 , (69) 



at the 99% CL [113]. It is clear that for n = 2 existing experiments are probing Mqg all the 
way up to the Planck scale M P1 ~ 10 19 GeV. For larger n, Eq.([68j) generalizes to 

HM" %o - 9xi(r ^ (70) 



This can be seen from Eq. (jobj) : the analysis of data places bounds on 7, and the neutrino 
energies are well above 1 GeV. Thus, for n — 3, the lower bound on Mqg > 10 13 GeV is 



34 



comparable to the one obtained from limits on proton decay 101]. Although the cubic 
energy dependence 7 oc E^ is not obtainable from the simple operator analysis presented 
above, it may be heuristically supported by a general argument that each of the bj must be 
suppressed by at least one power of M P1 , giving a leading behavior 

7 = k El/Ml, . (71) 

IceCube will collect a data set of order one million atmospheric neutrinos over 10 years. 
Not surprisingly, because of the increased energy and statistics over present experiments, 
the telescope will be sensitive to effects of quantum decoherence at a level well below current 
limits. Additionally, since the loss of quantum coherence is weighted by the distance traveled 
by the neutrinos, IceCube data analysis of high energy extraterrestrial neutrinos can be used 
to probe decoherence effects arising only after long propagation lengths. Such an analysis 
can be carried out by measuring the ratios of neutrino flavors present in the cosmic spectrum. 

Let the ratios of neutrino flavors at production in the cosmic sources be written as 
w e : Wfj, : w T with J2 a w a — 1, so that the relative fluxes of each mass eigenstate are given by 
Wj = ^2 a w a U^j- For ^> 1, the phases in Eq. (JoT]) will be erased by uncertainties in L 
and E v . Consequentl y, av eraging over sin 2 Ay one finds the transition probability between 
flavor states a and (3 115l 

P{v a ^v p ) = Y J Ul i U} i . (72) 

i 

Thus, we conclude that the probability of measuring on Earth a flavor a is given by 

P Va detected = ^ U *3 ^2 W ? U lj " ( 73 ) 

3 P 

Straightforward calculation shows that any initial flavor ratio that contains w e = 1/3 will 
arrive at Earth with equipartition on the three flavors. Since neutrinos from astrophysical 
sources are expected to arise dominantly from pions (and kaons) and their muon daughters, 
their initial flavor ratios of 1/3: 2/3:0 should arrive at Earth democratically distributed. 
Consequently there is a fairly robust prediction of 1/3 : 1/3 : 1/3 flavor ratios for measure- 
ment of astrophysical neutrinos. 

The prediction for the flavor population at Earth due to standard flavor-mixing (i.e. with 
no spacetime dynamics) of a pure z/ e beam is £\ l^plf/"^-! 2 ~ §|t^ a 2| 2 + §|^ai| 2 for flavor a, 
which leads to the flavor ratios ~ 3/5 : 1/5 : 1/5. This is very different from the democratic 
1/3 : 1/3 : 1/3. Since the effects of quantum decoherence would alter the flavor mixture to 
the ratio 1/3 : 1/3 : 1/3 (regardless of the initial flavor content) and the loss of quantum 
coherence is weighted by the distance travelled by the (anti) neutrinos, by measuring the 



Note that for k ~ 1, the lower limit on the UV cutoff is well beyond the electroweak scale. Thus, if 
TeV-scale gravity is realized in nature, interactions with virtual black holes would be non-dissipative. 
Non-dissipative interactions are expected when gravity is embedded in string theory, so that an S'-matrix 
description is possible. The existence of an S'-matrix makes it no longer automatic that, e.g., the B- 
violating ql and i?-conserving qq outgoing channels have the same probability, as they would in thermal 
evapora tion. Thus, the problem of avoiding rapid baryon decay is shifted to the examination of sym- 
metries |ll4) in the underlying string theory which would suppress the appropriate non-renormalizable 
operators at low energies. 
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i/-Cygnus beam IceCube will be sensitive to the effects of quantum decoherence at a level 



well below current limits 1161 . 



The Lindblad operators of an iV-level quantum mechanical system can be expanded in 
a basis of matrices satisfying standard commutation relations of Lie groups. For a 3-level 
system, the basis comprises the eight Gell-Mann SU(3) matrices plus the 3x3 identity 
matrix. As mentioned above, the theoretical approach provided by Lindblad quantum dy- 
namical semi-groups is a very general in the sense that no explicit hypothesis has been made 
about the actual interactions causing the loss of coherence. In what follows we adopt an 
expansion in a 3 flavor basis with a diagonal form for the 9x9 decoherence matrix, T>. Note 
that neutrinos oscillate among flavors separately between particle and antiparticle sectors 
and so the respective decoherence parameters for antineutrinos can be different from the 
corresponding ones in the neutrino sector. Upon averaging over the rapid oscillation for 
propagation between Cygnus OB2 and the Earth, only the diagonal Gel l-Ma nn matrices 



survive, and so the transition probability for antineutrinos takes the form [117 
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(74) 



where 7 3 and 7 8 are eigenvalues of the decoherence matrix for antineutrinos. Note that in 
Eq. ([71jl we set the CP violating phase to zero, so that all mixing matrix elements become 
real. Furthermore, under the assumptions that CPT is conserved and that decoherence 
effects are negligible at present experiments, the values of the mixing angle combinations 
appearing in Eq. (|74j) can be well determined by the usual oscillation analysis of solar, 
atmospheric, long-baseline and reactor data [8(1 ] . In what follows, we will assume that CPT 
is conserved both by neutrino masses and mixing as well as in decoherence effects. Note 
however that since the decoherence effects in the present study primarily affect antineutrinos, 
the result will hold for the antineutrino decoherence parameters exclusively if CPT is violated 
only through quantum-gravity effects. 

Now, we require further 7 3 = 7 8 = 7 ( = 73 = 7 8 under CPT conservation) so that 
Eq. (|74|) can be rewritten for the case of interest as: 
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(75) 



We make this simplification only to emphasize the primary signature of quantum decoher- 
ence, namely that after travelling a sufficiently long distance the flavor mixture is altered 
to the ratio 1/3 : 1/3 : 1/3, regardless of the initial flavor content. Consequently, if a 
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flux of antineutrinos were to be observed from the Cygnus spiral arm with a flavor ratio 
7^ 1/3 : 1/3 : 1/3, strong constraints can be placed on the decoherence parameter 7. 

Using the results of the up-to-date 3-v oscillation analysis of solar, atmospheric, long- 
baseline and reactor data |8(| we obtain the following values and 95% confidence ranges ^ 

Jve^vn — u - ±uu -0.082 ) 

Ur^r = -0.1281S , (76) 
t - n 0*7+0.011 

)v^v T — U.UJ/ _o.035 • 

The numbers given in Eq. ([7BJ) are obtained using the same techniques as described in 
Ref. but including the final SNO salt phase data |l!8j ]. 

Equipped with Eqs. (|75J). and (|75|) . we now proceed to determine the sensitivity of 

IceCube to decoherence effects. We estimate the expected number of induced tracks from 
the Cygnus OB2 source antineutrino flux using the semi-analytical calculation discussed 
in Sec. Ill Al Namely, we replace the antineutrino flux shown in Fig. |U] into Eq. (JBJ) and 
obtain, for standard mass- induced oscillations, a total of 212 x Pv e -*v — 48 V^- induced 
tracks in 15 years of observation (cluster within 1° of the source direction). For showers, 
the angular resolution is significantly worse than for muon tracks. Normally, a reduction of 
the muon bremsstrahlung background is effected by placing a cut of 4 x 10 4 GeV on the 
minimum reconstructed energy |8J. For Cygnus OB2, this strong energy cut is not needed 
since this background is filtered by the Earth. Thus we account for all events with shower 
energy E s ^ > E^ n = 10 3 GeV, trigger level. The directionality requirement, however, 
implies that the effective volume for detection of showers is reduced to the instrumented 
volume of the detector, 1 km 3 , because of the small size of the showers (less than 200 m in 
radius) in this energy range. Replacing the antineutrino flux shown in Fig. El into Eq. ()14j) . 
within the framework of standard oscillations, one expects 25 showers from the Cygnus 
OB2 source in 15 years of operation. We now turn to the estimate of the background. 
There are two different contributions - - atmospheric neutrinos and additional fluxes of 
extraterrestrial neutrinos. We obtain the number of expected track and shower events from 
atmospheric neutrinos as in Eqs ©, (|T5j). and ([TB^I with dF^" m /dE v as given in Fig. [HI being 
the u e and atmospheric neutrino fluxes integrated over a solid angle of of 1° x 1° (for 
tracks) and 10° x 10° (for showers) width around the direction of the Cygnus OB2 source 
[Q = 131.2°). We get an expected background of 14 atmospheric tracks and 47 atmospheric 
showers in 15 years. Of the 47 showers, 16 correspond to v e CC interactions while 31 
correspond to NC interactions. The large yield of NC events is due to the fact that at 
these energies, the atmospheric flux contains a very unequal mix of neutrino flavors (with 
ratios ~ 1/20 : 19/20 : 0). Interestingly, by increasing the minimum shower energy cut to 
5 x 10 3 GeV, v e CC and NC contribute in equal amounts to the number of at mosp heric 
showers. This is in agreement with simulations by the AMANDA Collaboration |l!9j . We 
turn now to the discussion of background events from other extraterrestrial sources. As 
discussed in Sec. 1111 CI the TeV 7-ray flux reported by the HEGRA Collaboration [H^ in the 
vicinity of Cygnus OB2 is likely to be due to hadronic processes. For the purposes of setting 
an upper bound on the neutrino flux we ignore all other sources near J2032+4130 because 
their steady emission in 7-rays is estimated to be smaller by more over a factor of 5 than the 
source of interest. 9 Using the expected flux of neutrinos from J2032+4130 given in Fig.|)]we 



This includes the famous microquasar Cygnus X-3, for which HEGRA Collaboration reported [^3 an 
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obtain the corresponding background from neutrinos with flavor ratios 1/3 : 1/3 : 1/3. As 
can be seen in Fig. the background is dominated by atmospheric neutrinos. Thus after 
15 years of data collection we expect 18 tracks and 1 shower from J2032+4130 for standard 
oscillations. 

We will now discuss how to isolate the possible signal due to decoherence in the antineu- 
trinos from Cygnus OB2 from the atmospheric background and possible fluctuations in the 
event rate due to unknown diffuse fluxes of extraterrestrial neutrinos. In general, we can 
predict that the expected number of track and shower events in the direction of the Cygnus 
OB2 source to be 

Aft? = AC + AC m + A" t HEGRA , (77) 
A4 h = Afl + AOT + Afl EGRA . (78) 

The first term corresponds to antineutrinos from neutron /3-decay. In the presence of deco- 
herence effects these event rates can be computed from Eqs. (JUJ), (|15|) and ()16|) with flavour 
transition probabilities given in Eq. (J75|) with L =1.7 kpc. The second term refers to atmo- 
spheric (anti)neutrinos (A/^ tm — 14, Af^ m = 47 for 15 years of exposure). The third term 
takes into account additional contributions from a diffuse flux of (anti)neutrinos produced 
via charged pion decay. In principle, decoherence effects may also affect the expected number 
of events from this diffuse flux. However given that the flavour ratios both from oscillation 
and decoherence are very close to 1/3 : 1/3 : 1/3 for the case of neutrinos produced via 
charged pion decay, there is no difference in the sensitivity region if decoherence effects are 
included or not in the evaluation of A/" t ™ GRA and A/" S REGRA . They are normalized to the 
maximum expected flux from J2032+4130 by a factor x = A/£ EGRA /18 = A/^ EGRA /1. 

Altogether, the quantities Atr and ACh, as defined in Eqs. (J77J) and (|7Sjl. can be regarded 
as the theoretical expectations of these event rates, corresponding to different points in the 
x — K n parameter space. For a given set of observed rates, A/" t ° bs and Af^ s , two curves 
are obtained in the two-dimensional parameter space by setting M°^ s = Af tr and Af°^ s = 
A/" s h- These curves intersect at a point, yielding the most probable values for the flux and 
decoherence scale for the given observations. Fluctuations about this point define contours 
of constant \ 2 i n an approximation to a multi-Poisson likelihood analysis. These contours 
can be obtained using Eq. (|4*B|) with the identification i = sh, tr. Marginalizing with respect 
to x, for n = 2 we obtain the following 1 degree-of-freedom bound 

7 < 2.0 x 1(T 44 (5.5 x 1(T 42 ) GeV , (79) 

at the 90% (99%) CL [9j. This corresponds to an improvement of abo ut 17 orders of 
magnitude over the best existing bounds from the SK and K2K data |l07j . Moreover, for 
n = 3, the 1-degree freedom bound is 

7 < 3.0 x 10~ 47 (2.9 x 10~ 45 ) GeV (80) 

at 90% (99%) CL @. For M QG = M P1 this corresponds to k < 3 x 10~ 3 at the 99% CL, well 
below the natural expectation, R ~ 1. 

Finally we note that, for 7 = k (E u /GeV)~ 1 , a non-vanishing 7 in Eq. (f?3|) can be related 
to a finite V e lifetime in the lab system |l20l | 

= g-i/nab — e ~Lm 17e /E 17 T 17e ^ 



upper limit for steady emission of F 1 (E 1 > 0.7 TeV) = 1.7 x 10 13 cm 2 s 1 . 
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where r Ve is the antineutrino rest frame lifetime and m Ve its mass. By duplicating our 
discussion for 7 oc E we obtain a 90% (99%) CL sensitivity 

7 < 1.0 x 10~ 34 (2.3 x 10~ 31 ) GeV , (82) 

which can be translated into 

> 10 34 GeV" 2 = 6.5 s eV" 1 . (83) 
This corresponds to an improvem ent o f about 4 orders of magnitude over the best existing 



bounds from solar neutrino data [121|, and of cours e giv es results comparable to the reach 
derived for neutrinos decaying over cosmic distances [122]. It should be noted that although 
the similar algebraic structure of the decoherence term in Eq. (|75|) and a decaying component 
in the neutrino beam provide a bound on the neutrino lifetime, these are conceptually two 
different processes. The decoherence case can be viewed as a successive absorption and re- 
emission of a neutrino from the quantum foam with change in flavor but no change in the 
average energy because of the condition [H, bj] = 0. This contrasts with the decay process 
which involves the emission of an additional particle. 



In summary, IceCube will provide a major improvement in the sensitivity to possible 
effects of quantum gravity. Because of the increase energy and statistics over present ex- 
periments measuring the atmospheric neutrino flux, existing bounds on the violation of the 
equivalence principle and of Lorentz invariance can be improved by two orders of magni- 
tude. On the other hand, antineutrinos produced in the decays of neutrons from Cygnus 
OB2 provide an excellent source in which to search for decoherence effects. Although the 
precise conclusions depend on the model considered, IceCube will improve the sensitivity 
to decoherence effects of 0{E 2 /Mpi) by 17 orders of magnitude over present limits and, 
moreover, it can probe decoherence effects of 0(E 3 /Mpj) which are well beyond the reach 
of other experiments. 



VI. DARK MATTER 



Over the past few years, a flood of h igh-q uality data from the Supernova Cosmology 
Project )l23| . the Supernova Search Team |l24j . the Wilkinson Microwave Anisotropy Probe 
fWMAP) jl25|. the Two Degree Field Galaxy Redshift Survey (2dFGRS) and the 

Sloan Digital Sky Survey (SDSS) |l27l | pin down cosmological parameters to-percent level 
precision, establishing a new paradigm of cosmology. A surprisingly good fit to the data 
is provided by a simple geometrically flat Friedman-Robertson- Walker Universe, in which 
30% of the energy density is in the form of non-relativistic matter (Q m = 0.30 ± 0.04) and 
70% in the form of a new, unknown dark energy component with strongly negative pressure 
(^a = 0.70 ± 0.04) jl28[. The matter budget has only 3 free parameters: the Hubble 
parameter h = 0.70^0'q 4 ,, the matter density Vt m h 2 = 0.138 ± 0.012, and the density in 
baryons, Vt^h 2 = 0.02301o;qq^- 10 This implies that the structure of the Universe is dictated 



The latter is consistent with the estimate from Big Bang nucleosy nthes is , based on measurements of 
deuterium in high redshift absorption systems, = 0.020 ± 0.002 |l29| . 
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FIG. 9: Simulations of structure formation with varying amounts of matter in the neutrino 
component, i.e., varying neutrino mass: (top left) massless neutrinos, (top right) ^2 a m Ua = 1 eV, 
(bottom left) ^2 a m Ua = 7 eV, (bottom right) Y1 Q = 4 eV. 



by the physics of as- yet-undiscovered cold dark matter (^cdm^ 2 = 0.115 ± 0.012) and the 
galaxies we see today are the remnants of relatively small overdensities in the nearly uniform 
distribution of matter in the very early Universe. Overdensity means overpressure that drives 
an acoustic wave into the other components making up the Universe: the hot gas of nuclei 
and photons and the neutrinos. These acoustic waves are seen today in the temperature 
fluctuations of the microwave background as well as in the distribution of galaxies on the 
sky. 

The only dark matter particle which is known to exist from experiment is the neutrino. 
Neutrinos decouple from thermal equilibrium while still relativistic, and consequently con- 
stitute the hot dark matter. With a contribution to the Universe's matter balance similar 
to that of light, neutrinos play a secondary role. The role is however identifiable: be- 
cause of their large mean free path, neutrinos prevent the smaller structures in the cold 
dark matter from fully developing and this is visible in the observed distribution of galax- 
ies, see Fig. [§] Simulations of structure formation with varying amounts of matter in the 
neutrino component can match to a variety of observations of today's sky, including mea- 
surement of galaxy-galaxy correlations and temperature fluctuations on the surface of last 
scattering. The resu lts suggest a neutrino mass, summed over the three neutrino flavors, 
^ a m„ o < 1 eV |l30j ]. a range compatible with the one deduced from oscillations. 

The si mple st model for cold dark matter consists of WIMP's - weakly interacting massive 



particles [1311 ]. Generic WIMP's were once in thermal equilibrium, but decoupled while 
strongly non-relativistic. Many approaches have been developed to attempt to detect dark 
matter. Such endeavors include direct detection experiments which hope to observe the 
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scattering of dark matter particles with the target material of the detector |132j and indirect 
detection experiments which are designed to s earch for the products of WIMP annihilation 
into gamma-rays, anti-matter and neutrinos 133l |. Neutrino telescopes indirectly search 



for the presence of dark matter by taking advantage of the Sun's ability to capture large 
numbers of WIMP's over time. Over billions of years, a sufficiently large number of WIMP's 
can accumulate in the Sun's core to allow for their efficient annihilation. Such annihilations 
produce a wide range of particles, most of which are quickly absorbed into the solar medium. 
Neutrinos, on the other hand, may escape the Sun and be detected in experiments on the 
Earth. The prospects for such experiments detecting dark matter critically depend on the 
capture rate of WIMP's in the Sun, which in turn depends on the elastic scattering cross 
section of these particles. 

The rate at which WIMP's are captured in the Sun depends on the nature of the inter- 
action the WIMP undergoes with nucleons in the Sun. These elastic scattering processes 
are often broken into two classifications: spin dependent interactions in which cross sec- 
tions increase with the spin of the target nuclei, and spin independent interactions in which 
cross sections increase with the total number of nucleons in the target. For spin-dependent 
interactions, the WIMP capture rate in the Sun is given by 1341 



where 



Cg^-S.SSxlO^S-V^local^locahmwiPM) Qpfj^) > ( 84 ) 
,, _ s ( Plocal \ / ^local \ ^ ( "iWIMP \ ~ 2 , OK x 

f (Piocai, tw, ™wipm) - ^ 3GeV/cm3 ) ^^m^ ) I lOOGeV ) ' (85) 

Piocai is the local dark matter density, <th,sd is the spin-dependent WIMP-on-proton (hy- 
drogen) elastic scattering cross section, viocai is the local rms velocity of halo dark matter 

particles, and m WIM p is the dark matter particle's mass. The analogous formula for the 

capture rate from spin-independent (scalar) scattering is jl34j 



n& o oc in20 -1 ft - \ ( °"H,SI + 0.07er H c,si\ , OCi , 

~ 3.35 x 10 s /(piocai, uiocai, ™wipm) I 1Q6 I, (86) 

where cth,si is the spin-independent WIMP-on-proton elastic scattering cross section and 
c"He,si is the spin-independent WIMP-on-helium, elastic scattering cross section. Typically, 
c"He,si — 16.0cth,si- The factor of 0.07 reflects the solar abundances of helium relative to 
hydrogen as well as dynamical factors and form factor suppression. Note that these capture 
rates are suppressed by two factors of the WIMP mass. One of these is simply the result of 
the depleted number density of WIMP's (n oc fn^ 1MP ) while the second factor is the result 
of kinematic suppression for the capture of a WIMP much heavier than the target nuclei, in 
this case hydrogen or helium. 

If the capture rate and annihilation cross section is sufficiently high, equilibrium will be 
reached between these two processes. For N WIMP's in the Sun, the rate of change of this 
quantity is given by 

^L = C e -A 6 N\ (87) 
at 

where C© = C® D +C® and A & = (av)/V cS . Here, V cS = 5.7x 10 27 cm 3 (m W iMp/100 GeV)" 3 / 2 , 
is the effective volume of the core of the Sun, which is estimated by matching the core tem- 
perature with the gravitational potential energy of a single WIMP at the core radius [l35j. 
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The present WIMP annihilation rate is given by 



^A & N 2 = ^C Q tanh 2 (VC^A^ * ) , (88) 



where t Q ~ 4.5 billion years is the age of the solar system. The annihilation rate is maximized 
when it reaches equilibrium with the capture rate. This occurs when 

VC & A®t Q > 1 . (89) 

For many particle dark matter candidates, this condition can be met. If this is the case, 
the final annihilation rate (and corresponding event rate or neutrino flux) has no further 
dependence on the dark matter particle's annihilation cross section. 

The sensitivity of direct detection experiments has been improving at a steady rate. The 
Cold Dark Matter Search (CDMS) experiment, operating in the Soudan mine in northern 
Minnesota, currently has produce d the strongest limits on spin-independent scattering cross 
sections of WIMP's with nucleons |l36j ]. as well as on spin-dependent scattering cross sections 
of WIMP's with neutrons |l37j |. CDMS data exclude spin-independent cross sections larger 
than approximately 2 x 10" 7 pb for a 50-100 GeV WIMP or 7 x 10" 7 pb (m W i M p/500 GeV) 
for a heavier WIMP 136]. The Zeplin-I |l38l | and Edelweiss |l39j experiments currently have 



spin-independent bo unds that are roughly a factor of 5 weaker over this mass range. The 
NAIAD experiment 140l | has placed the strongest constraints on spin-dependent WIMP- 



proton scattering. The data limit the spin-dependent cross section with protons to be 
less than approximately 0.3 pb for a WIMP in the mass range of 50-1 00 GeV and less 
than 0.8 pb (m WIMP /500 GeV) for a heavier WIMP. The PICASSO [l4l| and CDMS [HJ 
experiments have placed limits on the spin- dependent WIMP-proton cross section roughly 
one order of magnitude weaker than the NAIAD bound. All of these results (of course) 
impact t he pr ospects for detecting neutrinos produced by the annihilation of WIMP's in 
the Sun 142l |. Currently, the SK experiment has placed the strongest bounds on high- 



energy neutrinos from the direction of the Sun |143| . SK has two primary advantages over 
other experiments. Firstly, they have analyzed data over a longer period than most of their 
competitors, a total of nearly 1700 live days. Secondly, SK was designed to be sensitive 
to low energy (~ GeV) neutrinos, which gives them an advantage in searching for lighter 
WIMP's. SK's limit on neutrino-induced muons above 1 GeV from WIMP annihilations in 
the Sun is approximately 1000 to 2000 per square kilometer per year for WIMP's heavier 
than 100 GeV, and approximately 2000 to 5000 per square kilometer per year for WIMP's 
in the 20 to 100 GeV range. The precise va lue o f these limits dep ends on the WIMP 
annihilation modes considered. The AMANDA 144l | and MACRO 145l | collaborations have 



reported limits on the flux of neutrino-induced muons from the Sun that are only slightly 
weaker than SK. The limit placed by the AMANDA experiment resulted from only 144 live 
days of data. Having operated the detector for five years, AMANDA is expected to produce 
significantly improved bounds in the future. 

The supersymmetric extension of the Standard Mod el is the leading candidate to avoid 
't Hooft naturalness problem with the Higgs mass 1461 ] . Supersymmetry (SUSY) posits a 



"complete democracy" between integral and half-integral spins, implying the existence of 
many as- yet-undiscovered superpartners. Thus, if SUSY can serve as a theory of low energy 
interactions, it must be a broken symmetry. The most common assumption is that the 
minimal low energy effective supersymmetric theory (MSSM) has a breaking scale of order 
Asusy ~ 1 TeV. MSSM has a concrete advantage in embedding the Standard Model in a 
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GUT: the supersymmetric beta functions for extrapolating the measured strengths of the 
strong, electromagnetic, and wea k cou plings lead to convergence at a unified energy value 
of the order M GU t ~ 10 16 GeV |l47| . Moreover, by imposing i?-parity conservation one 
obtains as a byproduct the sta bility of the lightest SUSY particle x, making it a possible 
candidate for cold dark matter |l48j |. 

T he re lic abundance of SUSY WIMP's can be found by integrating the Boltzman equa- 
tion |l4fl| . 

d ^ + 3Hn = -(av)(n*-nl q ) J (90) 

where n is the present number density of SUSY WIMP's, n eq = g (m^ T/2n) 3 / 2 e~ m ^l T is 
the equilibrium number density, H is the expansion rate of the Universe at temperature T, 
g is the number of internal degrees of freedom of the WIMP, (av) is the thermally averaged 
annihilation cross section, and is the neutralino mass. Note that in the very early 
Universe, when n ~ n eq , the right hand side of Eq. (j9~Uj) is small and the evolution of the 
density is dominated by Hubble expansion. As the temperature falls below m%, however, the 
equilibrium number density becomes suppressed and the annihilation rate increases, rapidly 
reducing the number density. Finally, when the number density falls enough, the rate of 
depletion due to expansion becomes greater than the annihilation rate, H > n(av), and 
the neutralinos freeze out of thermal equilibrium. The freeze-out temperature Tf depends 
logarithmically upon (av), but for models with TeV scale SUSY breaking, one finds that 
Tf/m^ ~ 0.05. SUSY WIMPs are, by far, the favored candidate for cold dark matter, 
because for mas ses o f order 100 GeV to 10 TeV a present density of Q% h 2 ~ 0.1, comes out 
fairly naturally |l50j . n In what follows we mantain our discussion as general as possible but 
keep in mind the specific needs of the MSSM. 

As they annihilate, WIMPs can generate neutrinos through a wide range of channels. 
Annihilations to heavy quarks, tau leptons, gauge bosons and Higgs bosons can all generate 
neutrinos in the subsequent decay |l52j]. The total flux of neutrinos emitted by the Sun due 
to WIMP annihilition is then 



dE v 



dE, 



e -£Vl50 GeV ; (Q1) 



v , j 



where F^ q is the non-equilibrium suppression factor (« 1 for capture-annihilation equilib- 
rium ) , (d N„ / dE v )j is the flux produced by the j decay channel per WIMP annihila- 
tion 152l | . No te that neutrinos produced near the center of the Sun interact with the solar 
medium |l53| . yielding a depletion of the emission spectrum by a factor ~ e -E„/i5o Gev_ 
Finally, the flux reaching the Earth is, 

^ = ^-w- • (92) 

«±/(U LLlu v „„., ri , p 



An alternative dark matter candidate, which has received quite some attention recently, arises in models 
where all of the particles of the Standard Models can propagate through the bulk of an extra dimensional 
space that is compactified on a scale around ~ TeV -1 . Within this set up each Standard Model particle 
is accompanied by a tower of Kaluza Klein states. The lightest KK particle can be naturally stabl e in a 
way anologous to how the lightest superpartner is stabilized in i?-parity conserving models of SUSY |l5lj . 
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FIG. 10: IceCube event rate from dark matter annihilations in the Sun, as a function of the 
WIMP's spin-dependent elastic scattering cross section. In the left frame, no points shown violate 
the current spin-independent scattering constraints of C DMS. In the right frame, no points would 
violate the a spin- independent bound 100 times stronger [15J]. 



where d ~ 1.5 x 10 8 km is the Earth-Sun distance. 

By replacing this flux in Eq. (jUJ) it is easily seen that the CDMS bound implies that the 
event rate at IceCube for spin-independent interactions is at the level of the atmospheric 
neutrino background 1541 . Thus, isolating a signal from WIMP's which scatter with nu- 



cleons mostly spin- independently in the IceCube data sample would be a challenge. On 
the other hand, a WIMP with a largely spin-dependent scattering cross section with pro- 
tons may be capable of generating large event rates at IceCube. For example, a 300 GeV 
WIM P w ith a cross section near the experimental limit, leads to rates as high as ~ 10 6 per 
year |l54j . The elastic scattering and annihilation cross sections of a neutralino depend on 
its various couplings and on the mass spectrum of the Higgs bosons and superpartners. The 
neutralino couplings depend on its composition. Generally, the lightest neutralino can be 
any mixture of bino, wino and the two CP-even higgsinos, although in most models a largely 
bino-like neutralino is lightest. Spin dependent, axial- vector, scattering of neutralinos with 
quarks within a nucleon is made possible through the t-channel exchange of a Z, or the s- 
channel exchange of a squark. Spin independent scattering occurs at the tree level through 
s-channel squark exchange and t-channel Higgs exchange, and at the one-loop level through 
diagrams involving a loop of quarks and/or squarks. For generality, here we do not assume 
any particular SUSY breaking scenario or unification scheme, and scan the entire MSSM 
parameter space using the DarkSUSY program |l55| . We find that, after appling the cuts 
to comply with all collid er co nstraints and to produce a thermal relic density that saturates 
the observational limits |l56| | . very large spin-dependent cross sections (o"sd > 10 -3 pb) are 
possible even in models with very small spin-independent scattering rates. 

In Fig. El we show the event rate at IceCube from WIMP annihilation in the Sun versus 
the WIMP's spin-dependent cross section with protons. In the left frame, all points evade 
current constraints from CDMS. In the right frame we show the same result, but only showing 
those points which evade a constraint 100 times stronger than current CDMS bound. We 
can thus conclude that the next generation of direct detection experiments will not be able 
to probe the MSSM parameter space accessible to IceCube. 
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In summary, direct detection dark matter experiments, lead by the CDMS Collabora- 
tion, have placed increasingly stronger constraints on the cross section for elastic scattering 
of WIMP's on nucleons. Although these constraints impact the prospects for indirect de- 
tection experiemnts, the latest bounds placed on the scalar WIMP-nucleon cross section by 
the CDMS experiment do not seriously limit the WIMP parameter space to be probed by 
IceCube. This is largely due to the fact that the spin- dependent scattering of WIMP's with 
protons is the most efficient process for capture in the Sun for many particle dark matter 
models. The spin-dependent scattering cross section of a WIMP is not nearly as strongly 
constrained as the spin-independent quantity and large values of spin-dependent cross sec- 
tions are beyond the reach of any planned direct detection experiments. Thus, IceCube will 
play an important as well as complementary role in the search for particle dark matter. 



VII. TOPOLOGICAL DEFECTS AND SUPERHEAVY RELICS 

A plethora of explanations have been proposed to address the production mechanism of 
ultra-high energy cosmic rays |l57j |. In the absence of a single model which is consistent 
with all data, the origin of these particles remains a mystery. Clues to solve the mystery are 
not immediately forthcoming from the data, particularly since various experiments report 
mutually inconsistent results. 

In recent years, a somewhat confused picture vis-a-vis the energy spectrum and arrival 
direction distributi on h as been emerging. Since 1998, the AGASA Collaboration has con- 
sistently reported 158l | a continuation of th e spe ctrum beyond the GZK cutoff 23]. In 



contrast, the most recent results from HiRes [l59j describe a spectrum which is consistent 
with the expected GZK feature. This situation exposes the challenge posed by sy stem atic 
errors (predominantly arising from uncertainties in hadronic interaction models |l60l |) in 
these types of measurements. Although there seems to be a remarkable agreement among 
experiments on predictions about isotropy on large scale structure, this is certainly not 
the case when considering the two-point correlation function on a small angular scale. The 
AGASA Collaboration reports observations of event c luste rs which have a chance probability 



smaller than 1% to arise from a random distribution .161]. Far from confirming the AGASA 
results, the recent analysis reported by the HiRes Collaboratio n sho wed that their data are 
consistent with no clustering among the highest energy events 162j |. The discovery of such 



clusters would be a tremendous breakthrough for the field, but the case for them is not yet 
proved. Special care must be taken when computing the statistical significance in such an 
analysis. In particular, it is important to define the search procedure a priori in order to en- 
sure one does not inadvertently perform "trials" by studying the data before deciding upon 
the cuts. Very recently, with the aim of avoiding accidental bias on the number of trials 



performed in selecting the angular bin, the original claim of the AGASA Collabo ratio n 161 



was re-examined considering only those events observed after the original claim 163l |. This 



study showed that the evidence for clustering in the AGASA data set is weaker than was 
previously supposed, and is consistent with the hypothesis of isotropically distributed arrival 
directions. 

Further confusing the issue, recent HiRes data have been interprete d as a change in 



cosmic ray composition, from heavy nuclei to protons, at ~ 10 9 GeV 1641 ] . This is an 
order of m agnitude lower in energy than the previous crossover deduced from the Fly's Eye 
data |l65| . The end-point of the galactic flux is expected to be dominated by iron, as the 
large charge Ze of heavy nuclei reduces their Larmor radius (containment scales linearly with 
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Z) and facilitates their acceleration to highest energy (again scaling linearly with Z). The 
dominance of nuclei in the high energy region of the Galactic flux carries the implication 
that any changeover to protons represents the onset of dominance by an extra-galactic 
component. The inference from this new HiRes data is therefore that the extra-galactic flux 
is beginning to dominate the Galactic flux already at ~ 10 9 GeV. Significantly, this is well 
below E G zk ~ 10 10 7 GeV |23j], and so samples sources even at large redshift. 

The Pierre Auger Observatory is confronting the low statistics problem at the highest 
energies by instrumenting a huge collection area covering 3000 square kilometers on an ele- 
vated plane in Western Argentina The instrumentation consists of 1600 water Cerenkov 
detectors spaced 1.5 km apart. Showers occurring on clear moonless nights (about 10% 
of the operational time) are also viewed by four fluorescence detectors, allowing powerful 
reconstruction and cross-calibration techniques. Simultaneous observations of showers using 
two distinct detector methods will help to control the systematic errors that have plagued 
cosmic ray experiments to date. While no breakthroughs were reported |l66| . preliminary 
data forecast rapid progress and imminent results in deciding event rates at the high end of 
the spectrum. 

The difficulties so far encountered in modeling the production of ultra-high energy cosmic 
rays arise from the need to identify a source capable of launching particles to extreme en- 
ergy |l67j |. In contrast to the "bottom-up" acceleration of charged particles, the "top-down" 
scenario avoids the acceleration problem by assuming that charged and neutral pri marie s 
simply arise in the quantum mechanical decay of supermassive elementary X particles 168]. 
To maintain an appreciable decay rate today, it is necessary to tune the X lifetime to be 
longer (but not too much longer) than the age of the universe, or else "store" short-lived X 
particles in topological vestiges of early universe phase transitions. 

According to current unified models of high energy interactions, the Universe may have 
experienced several spontaneous symmetry breakings, where some scalar field, generally re- 
ferred to as the Higgs field, acquired a non-vanishing expectation value in the new vacuum 
(ground) state. Quanta associated with these fields are typically of the order of the GUT 
symmetry-breaking scale. During a phase transition, non-causal regions may evolve towards 
different states, so that in different domain borders the Higgs field may keep a null expec- 
tation value. Energy is then stored in a topological defect whose characteristics dep end on 
the topology of the manifold where the Higgs potential rea ches its minimum |l69^ . The 
relic de fects such as magnetic monopoles [170], cosmic strings 171 1, superconducting cosmic 
strings [1.7_2|, and cosmic necklaces [l73j (a possible hybr id to pological defect consisting of a 
closed loop of cosmic string with monopole "beads" on it H3) are all relatively topologically 
stable, but can release part of their energy (through radiation, annihilation, or collapse) in 
the form of X particles that typically decay to quarks and leptons. 

The highest energy cosmic rays may also be produced from the decay of some metastable 
superhea vy re l ic pa rticle with mass mx > 10 12 GeV and lifetime exceed ing the age of the 
Universe EH ElQ . Discrete gauged symmetries |l77j ] or hidden sectors |l7S| are generally 
introduced to stabilize the X particles. Higher dimensional operators, wormholes, and 
instanto ns ar e then inv o ked to break the new symmetry super-softly to maintain th e lon g 



lifetime |l75l . Il76l Il77l Il78j (collisional annihilation has been considered as well [179]) 



Arguably, these m etas t able super-heavy relics may constitute (a fraction of) the dark matter 
in galactic haloes 175l . 180j . If this were the case, due to the non-central position of the Sun 



in our Galaxy the fl ux of ultra-high energy cosmic rays from X particle decay would exhibit 
a dipole anisotropy |l8l| . 
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The cascade decay to cosmic ray particles is driven by the ratio of the volume density 
of the X-particle (% = p c Qx /m>x) to its decay time (t x ). This is very model dependent, 
as neither the cosmic average mass density contributed by the relics {&x), nor t x is known 
with any degree of confidence (p c « 1.05 x 10 _4 /i 2 GeV cm -3 ). Moreover, the internal 
mechanisms of the decay and the detailed dynamics of the first secondaries do depend on 
the exact nature of the particles. Consequently, no firm prediction on the expected flux of 
cosmic rays can be made. However, if there are no new mass scales between Agusy ~ 1 TeV 
and mx, 12 the squark and sleptons would behave like their corresponding supersymmetric 
partners, enabling one to infer from the "known" evolution of quarks and leptons the gross 
features of the X-particle cascade: the quarks hadronize pr oduc i ng je ts of hadrons containing 
mainly pions together with a 3% admixture of nucleons |l80l Il82j |. This implies that the 
injection spectrum is a rather hard fragmentation-type shape (with an upper limit usually 
fixed by the GUT scale) and dominated by 7-rays and neutrinos produced via pion decay. 
Therefore, the photon/proton ratio 183] and the ult ra-hi gh energy neutrino flux (which 
would be beyond the WB-flux of transparent sources) |l84j can be used as diagnostic tools 
to probe top-down models. 

In light of the mounting evidence that ultra- high energy cosmic rays are not 7-rays |l85| , 
one may try to force a proton dominance at ultra-high energies by postulating efficient 
absorption of the dominant ultra-high energy photon flux on the universal and/or galactic 
radio background. According to current estimates of the strengths of the magnetic field and 
of the radio wave background in our own galaxy, most ultra-high energy photons produced 
in the halo are expected to reach the Earth. In what follows we a ssum e that the photon 
interaction length in the galaxy has been greatly over-estimated jl86j |. and explore this 
assumption for neutrino signals. 13 

If a top down scenario is to explain the origin of ultra-high energy cosmic rays, the injec- 
tion spectrum should be normalized to account for the super-GZK events witho ut violating 
any observational flux measurements or limits at higher or lower energies |l88| . In partic- 
ular, neutrino and 7-ray fluxes depend on the energy released integrated over redshift, and 
thus on the specific top down model. Note that the electromagnetic energy injected into 
the Universe above the pair production threshold on the cosmic microwave background is 
recycled into a generic cascade spectrum below this threshold on a short time scale compared 
with the Hubble time. Therefore, it can have several potential observable effects, such as 
modified light element abundances due to 4 He ph otodisintegration, or induce spectral distor- 
tions of the cosmic microwave background |l89| . Additionally, measurements of the diffuse 
GeV 7-ray flux 38], to which the generic cascade spectrum would contribute directly, limit 
significantly the parameter space in whi ch cosmologically distant X's can generate the flux 
of the ultra-high energy cosmic rays 190j | . Such bounds are thus rel evan t only for topological 
defects, not for superheavy dark matter which is clustered locally |l86| | . 

A point worth noting at this juncture: as we discussed in the previous section, the cold 
dark matter hunt traditionally concentrates on WIMP's that were once in thermal equilib- 



12 This hypothesis, known as the desert "hypothesis" , is well motivated by the fact that the existence of new 
physics between the GUT scale and the SUSY breaking scale would destroy the very impresive feature of 
"natural" unification of the gauge couplings at Mgut ~ 10 16 GeV occuring in the MSSM. 

13 It is also possible that the energy measurements in air shower arrays are biased due to increasingly photon 
dominance at super- GZK energies, such that the superheavy dark matter model remains consistent with 
observational limits 11871. 
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FIG. 11: Right: The ultra-high energy cosmic ray flux predicted for the decay of superheavy 
particles with mass rax = 2 x 10 12 GeV is compared to the HiReS and AGASA cosmic ray 
data. The distribution of jets used includes an overdensity factor of 10 5 within 20 kpc of the 
galaxy. Spectra are shown for quark+antiquark (solid), quark+squark (dot-dash) and 5 quarks + 
5 squarks (dashes) initial states. Dark lines are from top-down origin alone whereas lighter lines are 
top-down plus contributions from an homogeneous population of astrophysical sources. Left: The 
neutrino plus antineutrino flux corresponding to the cosmic ray spectra from decay of superheavy 
particles 1921 ] . 



rium in the early Universe. Its present abundance is determined by the self-annihilation 
cross section. Since the largest annihilation cross section at early times is expected to be 
oc t^wimp; heavy WIMP's would have such a small cross section that their present abun- 
dance would be too large. Consequently, the mass of a thermal WIMP is found to be less 
than about 500 TeV. Then, in associating metastable superheavy relics with cold dark mat- 
ter we assume that the X-particles have never experienced thermal equil ibrium. Moreover, 
since the density of thermal WIMP's saturates the WMAP limit |l5q . the contribution 
from superheavy relics to the cold dark matter density must be much less than that of 
the relic neutralinos. As noted above the cosmic ray flux resulting from X-particle decays 
depends on the dimensionless parameter rx = (,xto/r x , where £x = ^x/^cdm, and to is 
the age of the universe. Scenarios which include in the rate normali zatio n photon flux from 
X-particles clustered in the halo lead to a value rx ~ 5 x 10~ n [175]. Omission of the 
photon channel in the normalization increases rx by a factor of about 2 — 10. Since models 
of X-production and decay typically lead to exponential dependence of both £ x and T x on 
a reheating temperature Tr and a quantum mechanical tunneling action, respectively; there 
is no impediment on accommodating this change in rx while maintaining Qx ^cdm- For 
example, assumuing a relic overdensity of 10 5 in the vici nity o f our galaxy, as expected for X 
particles that move freely under the influence of g ravit y [175j , an Earthly proton dominated 



spectrum requires £x oc e~ 2mx ^ TR ~ 10~ 4 — 10~ 8 191 . 

In Fig. ^2 we show the neutrino spectra produced in top-down models for a variety of X- 
decay channels. We have assumed that X particles have an overdensity of 10 5 in the vicinity 
of our galaxy. The exact profile of the halo of X particles does not affect the results as 
long as most ultra-high energy cosmic ray events originate at distances well below one GZK 
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X-decay channel 


AMANDA-II 


IceCube 




0.39 yr" 1 


12.2 yr- 1 


qq 


0.36 yr" 1 


11.4 yr" 1 


5 x qq 


1.40 yr" 1 


44.6 yr- 1 



TABLE II: Neutrino event rates at AMANDA-II and IceCube for various A-decay channels. 

interaction length. The spectra of stable particles (protons, photons, neutralino, elec tron s 
and neutrinos of the three species) have been computed using the SH-decay program [l93j. 
At the energies under consideration, it is necessary to take into account all the gauge cou- 
plings of the MSSM; indeed, at the scale of unification, they are all of the same strength, 
so that electroweak (and some Yukawa) interactions can be as relevant as the QCD ones. 
The perturbative part of the shower was computed by solving numerically the complete set 
of evolution equations for the relevant fragmentation functions of the MSSM. We modelled 
carefully the decays of unstable particles with mass near AgusY ~ 1 TeV, as well as the 
hadronization process at the GeV scale for light quarks and gluons. 14 The normalization is 
determined by matching the X-particle baryon flux to the difference between the observed 
cosmic ray spectrum at E ~ 10 11 GeV and contributions from a homogeneous population of 
astrophysical sources. All fluxes shown in Fig. ^2 are consistent with current measurements 
of the diffuse GeV 7-ray background. 

The event rates at AMANDA-II and IceCube are given in Table O [192]. They are 
estimated using Eqs. (jBJ) and ()14)) with a threshold energy 10 5 GeV for showers and muons. 
This energy cut effectively removes any background from atmospheric neutrino events. For 
simplicity, the muon effective area and shower effective volume of the AMANDA-II detector 
are taken to be uniform: A® s ~ 0.05 km 2 j5?| and V e fr ~ 0.008 km 3 (H, respectively. 

In summary, IceCube can test the viability of top-down models. If a signal is found, 
future high statistics experiments should be able to map out the neutrino spectrum, therby 
allowing us experimental access to physics at energy scales many order of magnitude beyond 
the scope of any conceivable particle collider on Earth. 

On a separate track, the cubic kilometer of ice provides a large detector area for direct 
search of fast magnetic monopoles. This topological defect appears in phase transitions 
which leave an unbroken U(l) symmetry group. A symmetry-breaking temperature T c ~ 
(H) at which the order parameter (H) turns on, leads to a monopole mass mx m ~ T c /a, 
where a is the fine structure constant at scale T c . To avoid violations of Standard Model 
physics the value of (H) should be at or above the electroweak scale, yielding a lower limit 
on the monopole mass: mx m > 40 TeV. 

The number density and therefore the flux o f monopoles emerging from a phase transition 
are determined by the Kibble mechanism [169]. At the time of the phase transition, roughly 
one monopole or antimonopole is produced per correlated volume. The resulting monopole 
number density today is 

n Xin ~ 10- 19 (T c /10 n GeV) 3 (W£ c ) 3 cm" 3 , (93) 
where £ c is the phase transition correlation length, bounded from above by the horizon 



The primary 10-body decay X — > bq5q has been modeled using phase space only, i.e. ignoring any possible 
dependence of the matrix element on external momenta. 
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size Ih at the time when the system relaxes to the true broken-symmetry vacuum. In a 
second order or weakly first order phase transition, the correlation length is comparable to 
the horizon size. In a strongly first order transition, the correlation length is considerably 
smaller than the horizon size. 

These monopoles easily pick up energy from the magnetic fields permeating the Universe 
and can traverse unscathed through the primeval radiation. The kinetic energy imp arted to 
a magnetic monopole on traversing a magnetic field along a particular path is ESI 



E = g B.dl ~g B £ y/E , (94) 

J path 

where g = e l '2a = 3.3 x 10~ 8 dynes/G is the magnetic charge according to the Dirac 
quantization [195J , B is the magnetic field strength, £ specifies the field's coherence length, 
and yfn is a factor to approximate the random-walk through the n domains of coherent 
fields traversed by the path. From Eq. ()93|) then, the general expression for the relativistic 
monopole flux may be written as 

4>x m = cn x j4n ~ 2 x 1(T 4 (^^Y (!)" cm" 2 sec" 1 sr" 1 . (95) 

Phenomenologically, the monopole flux is constrained by cosmology and by astrophysics. 
Cosmology requires that the monopole energy density Qx m should be less than the 
cold dark matter density of the Universe. From Eq. (|93|) comes a constraint ~ 
0.1 (m Xm /10 13 GeV) 4 (/ H /£ c ) 3 , which implies that monopoles with m Xm < 10 13 (£ C // H ) 3/4 GeV 
do not "over-close" the Universe. Astrophysics requires that there not be so many 
monopoles around a s to effectively "short out" the Galactic magnetic field, <px m < 
10~ 15 cm" 2 s _1 sr" - 1 Ell. To satisfy this constraint, the Kibble flux in Eq. (J93j) requires 



m x m < 10 11 (£ c /^) GeV . Note that with the usual GUT scale, the fractional monopole mass 
density fix m ~ 10 16 overdoses the Universe by sixteen orders of magnitude. However, to 



or 



dilute the monopole density inflation can be invoked after the phase transition 197 
else the U(l) group can be broken temporarily so as to create cosmic strin g defects which 
connect monopoles to anti-monopoles pairwise, which then annihilate |l98|. 

A magnetic monopole with velocity v — > c would emit Cerenkov radiation along its path. 
The total power emitted per unit frequency v and per unit length I 



dvdl 4a ' 1 } 



1 2 



v 2 n 2 



exceeds that of a bare relativistic muon by n 2 /4a, where n r is the refractive index of the 
medium 199l |. This corresponds to an enhancement factor of 4700 for monopoles interacting 



in vacuum and 8300 for monopole interactions in ice. Clear ly, th e large light output of a 
monopole track would be a rather unique signal at IceCube |200j. The non-observation of 
monopole signals at the large Cerenkov detector in the lake Baikal leads to a 90% CL upper 
limit = 0.5 - 0.7 - 1.9 x 10~ 16 cm^s^sr" 1 for p = 1.0 - 0.9 - 0.8 respectively 
(where p = v/c) HHI- The AMANDA (data taken during 1997) and MACRO p03j 
collaborations reported slightly weaker flux limits. Since the muon effective area of IceCube 
is about a factor of 20 larger than that of AMANDA, in 5 years of operation the kilometer 
scale telescope will reach a sensitivity nearly two orders of magnitude below existing limits. 
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VIII. OUTLOOK 



Today, precision data from man-made accelerators can, without exception, be accom- 
modated by the Standard Model of particle physics. Whenever the experimental precision 
increases, the higher precision measurements invariably collapse into Standard Model values. 
Nevertheless, since all the big questions remain unanswered, there is no feeling that we are 
now dotting the i's and crossing the t's of a mature theory. Worse, the theory has its own 
demise built into its radiative corrections. 

Led by a string of fundamental experimental measurements that unmasked the leptonic 
flavor mixing, neutrino physics has wounded the Standard Model. The hope is that at this 
point the glass is half full and that high-precision high-statistics IceCube data will continue 
the process and pierce the Standard Model's resistant armor. 
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